Cosmological parameters from SDSS and WMAP 
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We measure cosmological parameters using the three-dimensional power spectrum P(k) from 
over 200,000 galaxies in the Sloan Digital Sky Survey (SDSS) in combination with WMAP and 
other data. Our results are consistent with a "vanilla" flat adiabatic ACDM model without tilt 
(n s = 1), running tilt, tensor modes or massive neutrinos. Adding SDSS information more than 
halves the WMAP-only error bars on some parameters, tightening la constraints on the Hubble 
parameter from h m 0.74to'o7 to h m 0.70^0 03, on the matter density from fi m w 0.25 ± 0.10 to 
n m w 0.30 ± 0.04 (la) and on neutrino masses from < 11 eV to < 0.6 eV (95%). SDSS helps even 
more when dropping prior assumptions about curvature, neutrinos, tensor modes and the equation 
of state. Our results are in substantial agreement with the joint analysis of WMAP and the 2dF 
Galaxy Redshift Survey, which is an impressive consistency check with independent redshift survey 
data and analysis techniques. In this paper, we place particular emphasis on clarifying the physical 
origin of the constraints, i.e., what we do and do not know when using different data sets and prior 
assumptions. For instance, dropping the assumption that space is perfectly flat, the WMAP-only 
constraint on the measured age of the Universe tightens from to ~ 16.3 Gyr to to ~ 14.llo'g 
Gyr by adding SDSS and SN la data. Including tensors, running tilt, neutrino mass and equation 
of state in the list of free parameters, many constraints are still quite weak, but future cosmological 
measurements from SDSS and other sources should allow these to be substantially tightened. 



I. INTRODUCTION 

The spectacular recent cosmic microwave background 
(CMB) measurements from the Wilkinson Microwave 
Anisotropy Probe (WMAP) [1-7] and other experiments 
have opened a new chapter in cosmology However, as 



emphasized, e.g., in [6] and [8], measurements of CMB 
fluctuations by themselves do not constrain all cosmologi- 
cal parameters due to a variety of degeneracies in parame- 
ter space. These degeneracies can be removed, or at least 
mitigated, by applying a variety of priors or constraints 
on parameters, and combining the CMB data with other 



2 



cosmological measures, such as the galaxy power spec- 
trum. The WMAP analysis in particular made use of 
the power spectrum measured from the Two Degree Field 
Galaxy Redshift Survey (2dFGRS) [9-11]. 

The approach of the WMAP team [6, 7], was to ap- 
ply Ockham's razor, and ask what minimal model (i.e., 
with the smallest number of free parameters) is consis- 
tent with the data. In doing so, they used reasonable 
assumptions about theoretical priors and external data 
sets, which allowed them to obtain quite small error bars 
on cosmological parameters. The opposite approach is 
to treat all basic cosmological parameters as free param- 
eters and constrain them with data using minimal as- 
sumptions. The latter was done both in WMAP accu- 
racy forecasts based on information theory [12-16] and in 
many pre- WMAP analyses involving up to 11 cosmolog- 
ical parameters. This work showed that because of phys- 
ically well-understood parameter degeneracies, accurate 
constraints on most parameters could only be obtained 
by combining CMB measurements with something else. 
Bridle, Lahav, Ostriker and Steinhardt [8] argue that in 
some cases (notably involving the matter density fl m ), 
you get quite different answers depending on your choice 
of "something else" , implying that the small formal error 
bars must be taken with a grain of salt. For instance, 
the WMAP team [6] quote Q m = 0.27 ± 0.04 from com- 
bining WMAP with galaxy clustering from the 2dFGRS 
and assumptions about spatial flatness, negligible tensor 
modes and a reionization prior, whereas Bridle et al. [8] 
argue that combining WMAP with certain galaxy clus- 
ter measurements prefers tt m ~ 0.17. In other words, 
WMAP has placed the ball in the non-CMB court. Since 
non-CMB measurements are now less reliable and precise 
than the CMB, they have emerged as the limiting factor 
and weakest link in the quest for precision cosmology. 
Much of the near-term progress in cosmology will there- 
fore be driven by reductions in statistical and systematic 
uncertainties of non-CMB probes. 

The Sloan Digital Sky Survey [17-19] (SDSS) team has 
recently measured the three-dimensional power spectrum 
P(k) using over 200,000 galaxies. The goal of that mea- 
surement [20] was to produce the most reliable non-CMB 
data to date, in terms of small and well-controlled sys- 
tematic errors, and the purpose of the present paper is to 
use this measurement to constrain cosmological parame- 
ters. The SDSS power spectrum analysis is completely in- 
dependent of that of the 2dFGRS, and with greater com- 
pleteness, more uniform photometric calibration, analyt- 
ically computed window functions and improved treat- 
ment of non-linear redshift distortions, it should be less 
sensitive to potential systematic errors. We emphasize 
the specific ways in which large-scale structure data re- 
moves degeneracies in the WMAP-only analysis, and ex- 
plore in detail the effect of various priors that are put on 
the data. The WMAP analysis using the 2dFGRS data 
[6, 7] was carried out with various strong priors: 



1. reionization optical depth r < 0.3, 

2. vanishing tensor fluctuations and spatial curvature 
when constraining other parameters, 

3. that galaxy bias was known from the 2dFGRS bis- 
pectrum [21], and 

4. that galaxy redshift distortions were reliably mod- 
eled. 

We will explore the effect of dropping these assumptions, 
and will see that the first three make a dramatic differ- 
ence. Note in particular that both the spectral index n s 
and the tensor amplitude r are motivated as free parame- 
ters only by inflation theory, not by current observational 
data (which are consistent with n s = 1, r = 0), suggest- 
ing that one should either include or exclude them both. 

The basic observational and theoretical situation is 
summarized in Figure 1. Here we have used our Monte 
Carlo Markov Chains (MCMC, described in detail below) 
to show how uncertainty in cosmological parameters (Ta- 
ble 1) translates into uncertainty in the CMB and mat- 
ter power spectra. We see that the key reason why SDSS 
helps so much is that WMAP alone places only very weak 
constraints on the matter power spectrum P(k). As sim- 
plifying theoretical assumptions are added, the WMAP 
P(k) predictions are seen to tighten into a narrow band 
whose agreement with the SDSS measurements is a strik- 
ing manifestation of cosmological consistency. Yet even 
this band is still much wider than the SDSS error bars, 
which is why SDSS helps tighten constraints (notably on 
17a and h) even for this restricted 6-parameter class of 
models. 

The rest of this paper is organized as follows. After 
presenting our basic results in three tables, we devote a 
series of sections to digesting this information one piece 
at a time, focusing on what we have and have not learned 
about the underlying physics, and on how robust the var- 
ious conclusions are to the choice of data sets and prior 
assumptions. In Section VIII we discuss our conclusions 
and potential systematic uncertainties, assess the extent 
to which a robust and consistent cosmological picture 
emerges, and comment on upcoming prospects and chal- 
lenges. 
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FIG. 1: Summary of observations and cosmological models. Data points are for unpolarized CMB experiments combined (top; Appendix 
A. 3 details data used) cross-polarized CMB from WMAP (middle) and Galaxy power from SDSS (bottom). Shaded bands show the 1- 
sigma range of theoretical models from the Monte-Carlo Markov chains, both for cosmological parameters (right) and for the corresponding 
power spectra (left). From outside in, these bands correspond to WMAP with no priors, adding the prior /„ = 0, w = — 1, further adding 
the priors = r = a = 0, and further adding the SDSS information, respectively. These four bands essentially coincide in the top two 
panels, since the CMB constraints were included in the fits. Note that the £-axis in the upper two panels goes from logarithmic on the left 
to linear on the right, to show important features at both ends, whereas the fc-axis of the bottom panel is simply logarithmic. 
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II. BASIC RESULTS 
A. Cosmological parameters 

In this paper, we work within the context of a hot 
Big Bang cosmology with primordial fluctuations that 
are adiabatic (i.e., we do not allow isocurvature modes) 
and Gaussian, with negligible generation of fluctuations 
by cosmic strings, textures, or domain walls. Within this 
framework, we follow [6, 22] in parameterizing our cos- 
mological model in terms of 13 parameters: 

p = (T,u> b ,u) d ,fu,^A,w,fl k ,A s ,n s ,a,r,nt,b). (1) 

The meaning of these 13 parameters is described in Ta- 
ble 1, together with an additional 16 derived parameters, 
and their relationship to the original 13. 

All parameters are defined just as in version 4.3 of 
CMBfast [23]: in particular, the pivot point unchanged 
by n s , a and n t is at 0.05/Mpc, and the tensor normal- 
ization convention is such that r — — 8n t for slow-roll 
models, eg, the linear rms mass fluctuation in spheres of 
radius 8/i Mpc, is determined by the power spectrum, 
which is in turn determined by p via CMBfast. The 
last six parameters in the table are so-called normal pa- 
rameters [24] , which correspond to observable features in 
the CMB power spectrum [25, 26] and are useful for hav- 
ing simpler statistical properties than the underlying cos- 
mological parameters as discussed in Appendix A. Since 
current n t -constraints are too weak to be interesting, we 
make the slow-roll assumption n t = — r/8 throughout 
this paper rather than treat n t as a free parameter. 

B. Constraints 

We constrain theoretical models using the Monte Carlo 
Markov Chain method [27-33] implemented as described 
in Appendix A. Unless otherwise stated, we use the 
WMAP temperature and cross-polarization power spec- 
tra [1-4], evaluating likelihoods with the software pro- 
vided by the WMAP team [7]. When using SDSS in- 
formation, we fit the nonlinear theoretical power spec- 
trum P(k) approximation of [34] to the observations re- 
ported by the SDSS team [20], assuming an unknown 
scale-independent linear bias b to be marginalized over. 
This means that we use only the shape of the measured 
SDSS power spectrum, not its amplitude. We use only 
the measurements with k < 0.2ft,/Mpc as suggested by 
[20]. The WMAP team used this same /c-limit when an- 
alyzing the 2dFGRS [7]; we show in Section VIII C that 
cutting back to k < 0.15ft./Mpc causes a negligible change 
in our best-fit model. To be conservative, we do not 
use the SDSS measurement of rcdshift space distortion 
parameter [3 [20], nor do we use any other information 
("priors") whatsoever unless explicitly stated. When us- 
ing SN la information, we employ the 172 SN la redshifts 
and corrected magnitudes compiled and uniformly ana- 
lyzed by Tonry et al. [35], evaluating the likelihood with 



the software provided by their team, which marginalizes 
over the corrected SN la "standard candle" absolute mag- 
nitude. Note that this is an updated and expanded data 
set from that available to the WMAP team when they 
carried out their analysis [6]. 

Our constraints on individual cosmological parameters 
are given in Tables 2-4 and illustrated in Figure 2, both 
for WMAP alone and when including additional informa- 
tion such as that from the SDSS. To avoid losing sight 
of the forest for all the threes (and other digits), we will 
spend most of the remainder of this paper digesting this 
voluminous information one step at a time, focusing on 
what WMAP and SDSS do and don't tell us about the 
underlying physics. The one-dimensional constraints in 
the tables and Figure 2 fail to reveal important informa- 
tion hidden in parameter correlations and degeneracies, 
so a powerful tool will be studying the joint constraints 
on key 2-parameter pairs. We will begin with a simple 
6-parameter space of models, then gradually introduce 
additional parameters to quantify both how accurately 
we can measure them and to what extent they weaken 
the constraints on the other parameters. 



III. VANILLA ACDM MODELS 

In this section, we explore constraints on six-parameter 
"vanilla" models that have no spatial curvature (Qfc = 0) , 
no gravity waves (r — 0), no running tilt (a = 0), neg- 
ligible neutrino masses (/„ = 0) and dark energy cor- 
responding to a pure cosmological constant (w = — 1). 
These vanilla ACDM models are thus determined by 
merely six parameters: the matter budget (f^A, ^>d, ^b), 
the initial conditions (A s , n s ) and the reionization optical 
depth t. (When including SDSS information, we bring 
in the bias parameter b as well.) 

Our constraints on individual cosmological parameters 
are shown in Tables 2-4 and Figure 2 both for WMAP 
alone and when including SDSS information. Several fea- 
tures are noteworthy. 

First of all, as emphasized by the WMAP team [6] , er- 
ror bars have shrunk dramatically compared to the situ- 
ation before WMAP, and it is therefore quite impressive 
that any vanilla model is still able to fit both the un- 
polarized and polarized CMB data. The best fit model 
(Table 2) has x 2 ~ 1431.5 for 899 + 449 - 6 = 1342 ef- 
fective degrees of freedom, i.e., about 1.7er high if taken 
at face value. The WMAP team provide an extensive 
discussion of possible origins of this slight excess, and ar- 
gue that it comes mainly from three unexplained "blips" 
[7, 36] , deviations from the model fit over a narrow range 
of £, in the measured temperature power spectrum. They 
argue that these blips have nothing to do with features 
in any standard cosmological models, since adding the 
above-mentioned non- vanilla parameters does not reduce 
X 2 substantially — we confirm this below, and will not 
dwell further on these sharp features. Adding the 19 
SDSS data points increases the effective degrees of free- 
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Tabic 1: Cosmological parameters used. Parameters 14-28 arc determined by the first 13. Our Monte-Carlo Markov Chain assigns a 
uniform prior to the parameters labeled "MCMC" . The last six and those labeled "Fits" are closely related to observable power spectrum 
features [24-26] and are helpful for understanding the physical origin of the constraints. 



Parameter 


Meaning 


Status 


Use 


Definition 




T 


Reionization optical depth 


Not optional 










Baryon density 


Not optional 


MCMC 


u> b = Q b h 2 = p 6 /(1.88 X Kr 26 kg/m 3 ) 






Dark matter density 


Not optional 


MCMC 


ui d = U d h 2 = p d /(1.88 X Kr 26 kg/m 3 ) 




u 


Dark matter neutrino fraction 


Well motivated MCMC f v = p v / ' p d 




n A 


Dark energy density 


Not optional 


MCMC 






w 


Dark energy equation of state 


Worth testing 


MCMC pa/pa (approximated as constant) 




n k 


Spatial curvature 


Worth testing 








A s 


Scalar fluctuation amplitude 


Not optional 




Primordial scalar power at k = 0.05/Mpc 




n a 


Scalar spectral index 


Well motivated MCMC Primordial spectral index at k = 0.05/Mpc 




a 


Running of spectral index 


Worth testing 


MCMC 


a = d In n s jd In k (approximated as constant) 




r 


Tcnsor-to-scalar ratio 


Well motivated 


MCMC Tensor-to-scalar power ratio at k = 0.05/Mpc 




n t 


Tensor spectral index 


Well motivated 


MCMC 






b 


Galaxy bias factor 


Not optional 


MCMC 


b = [• p K alaxy( fc )/- p ( fc )] 1/2 (assumed constant for k < 0.2/i/Mpc) 


^ion 


Reionization redshift (abrupt) 






Zion ~ 92(0.03fer/wi,) 2 / 3 Slm 3 (assuming abrupt reionization; 


[37]) 




Physical matter density 




Fits 


UJm = Ul(> + Ul d = Q m h 2 






Matter density/critical density 






n m = i — Qa — f^fe 






Total density /critical density 






f2tot = n m + Qa = 1 — f2fc 




At 


Tensor fluctuation amplitude 






A t = rA s 




M v 


Sum of neutrino masses 






M v ai (94.4 eV) x ui d f v [38] 




h 


Hubble parameter 






h = \J (ojd + Mb)/ (i - ^fe - n A ) 




p 


Redshift distortion parameter 






P » + (1 + n m /2)(n A /70)]/6 [39, 40] 




to 


Age of Universe 






to ^(9.785 Gy^xfe- 1 f^[(Q A a-( 1 + 3 ^ +Q k + n m /a)]- 1 / 2 d 


i [38] 




Galaxy fluctuation amplitude 






a 8 = {4^/ °°[^(sinx - xcosx )] 2 P{k)^Y/\ x = kx8h 


_1 Mpc 


z 


CMB peak suppression factor 




MCMC 


Z = e' 2T 




Ap 


Amplitude on CMB peak scales 




MCMC 


Ap — As& ^ T 




e s 


Acoustic peak scale (degrees) 




MCMC e s (£lk,£lA,w,0Jd,0Jb) given by [25] 




H 2 


2nd to 1st CMB peak ratio 




Fits 


H 2 = (0.925^ 18 2.4 n =- 1 )/[l + K/0.0164) 12 "m' )]0- 2 [25] 




Hi 


3rd to 1st CMB peak ratio 




Fits 


#3 = 2.17[1 + (o; i ,/0.044) 2 ]- 1 ^ 59 3.6"=- 1 /[l + 1-63(1 - u, 6 /0.071)u; m ] 




Amplitude at pivot point 




Fits 


A, = 0.82 n =- 1 A p 





dom by 19 — 1 = 18 (since this requires the addition of 
the bias parameter b), yet raises the best- fit \ 2 by only 
15.7. Indeed, Figure 1 shows that even the model best 
fitting WMAP alone does a fine job at fitting the SDSS 
data with no further parameter tuning. 



A. The vanilla banana 

Second, our WMAP-only constraints are noticeably 
weaker than those reported by [6], mostly because we 
did not place a prior on the value of the reionization op- 
tical depth t, and adding SDSS information helps rather 
dramatically with all of our six basic parameters, roughly 
halving the 2a error bars. The physical explanation for 
both of these facts is that the allowed subset of our 6- 
dimcnsional parameter space forms a rather elongated 
banana-shaped region. In the 2-dimensional projections 
shown (Figures 3, 4, 5 and 6), this is most clearly seen 
in Figures 3 and 5. Moving along this degeneracy ba- 
nana, all six parameters (r, f^A, ^><i, ^b, A s , n s ) increase 



together, as does h. 

There is nothing physically profound about this one- 
dimensional degeneracy. Rather, it is present because 
we are fitting six parameters to only five basic observ- 
ables: the heights of the first three acoustic peaks, the 
large-scale normalization and the angular peak location. 
Within the vanilla model space, all models fitting these 
five observables will do a decent job at fitting the power 
spectra everywhere that WMAP is sensitive [25]. As 
measurements improve and include additional peaks, this 
approximate degeneracy will go away. 

Here is how the banana degeneracy works in practice: 
increasing r and A s in such a way that A p = A s e~ 2r 
stays constant, the peak heights remain unchanged and 
the only effect is to increase power on the largest scales. 
The large-scale power relative to the first peak can be 
brought back down to the observed value by increasing 
n s , after which the second peak can be brought back 
down by increasing u>b- Adding WMAP polarization in- 
formation actually lengthens rather than shortens the 
degeneracy banana, by stretching out the range of pre- 
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FIG. 2: Constraints on individual cosmological quantities using WMAP alone (shaded yellow/light grey distributions) and including SDSS 
information (narrower red/dark grey distributions). Each distribution shown has been marginalized over all other quantities in the class 
of 6-parameter (t, f2 a , > ^>b ; A a , n s ) "vanilla" models as well as over a galaxy bias parameter b for the SDSS case. The a-distributions 
are also marginalized over r and Oj.. The parameter measurements and error bars quoted in the tables correspond to the median and the 
central 68% of the distributions, indicated by three vertical lines for the WMAP+SDSS case above. When the distribution peaks near 
zero (like for r), we instead quote an upper limit at the 95th percentile (single vertical line). The horizontal dashed lines indicate e~ x / 2 
for x = 1 and 2, respectively, so if the distribution were Gaussian, its intersections with these lines would correspond to la- and 2a limits, 
respectively. 



ferred r- values — the largest-scale polarization measure- 
ment prefers very high r (Figure 1) while the unpolarized 
measurements prefer r — 0. This banana degeneracy was 
also discussed in numerous accuracy forecasting papers 



and older parameter constraint papers [12, 13, 15, 16]. 

Since the degeneracy involves all the parameters, es- 
sentially any extra piece of information will break it. 
The WMAP team break it by imposing a prior (assuming 
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Table 2: la constraints on cosmological parameters using WMAP information alone. The columns compare different theoretical priors 
indicated by numbers in italics. The penultimate column has only the six "vanilla" parameters (r, Q^, uj^, u>i, A s , n s ) free and therefore 
gives the smallest error bars. The last column uses WMAP temperature data alone, all others also include WMAP polarization 

information. 



Using WMAP temperature and polarization information 



6par+$7fc + r + a 6par+f2j; 



6par+r 



6par+/i, 



6par+«j 



6par 



No pol. 
6par 



e s 

n A 

h 2 n d 

h 2 n b 

U 

n s 

n t + 1 

Ap 

r 

b 

w 

a 



0.52 



+0.21 
0.15 



0.65 



-0.19 
-0.32 



0.68 



+0.13 
0.16 



0.75 



+0.12 
0.23 



0.68 



+ 0.15 
0.21 



0.66 



+0.17 
0.25 
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603+ ' 015 



5Q68+ ' 0048 

u.oyoo_ 0056 



=sq<j+0.0062 
u.ooyj_ 0056 
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u - 09 °' -0.0048 



0.54 



+0.24 
0.33 



0.105lg :u 23 3 
0.02381°;^? 





+0.13 
-0.10 
+0.0097 
0.0141 



0.53 
0.108 



0.97 
0.9847 

593+0 053 
u.ono_ 044 

< 0.90 (95%) 

No constraint 



+0.24 
0.32 
+0.022 
0.034 

0241+ - 0055 

U.UZ4tl_ 0020 





+ 0.18 
-0.06 

1 

+0.053 
0.051 



0.823 
0.097 



+0.058 
-0.082 
+0.021 
0.018 



0.687 



+0.087 
0.097 



0256+ ' 0025 

U.U^OD_q 00ig 





llq +0.018 
U.11M_q 016 

0247+ ' 0029 
u - uz *' -0.0016 



0.64 
0.118 



+0.14 
0.17 
+0.020 
0.020 



0246+ - 0038 
u.u^+o_ 0017 



+0.10 
0.10 
+0.020 
0.021 

0245+ - 0050 
u.uz+o_ 0019 



0.75 
0.115 



> 0.50 (95%) 

5q84 +0 ' 0041 
u - oys4 -0.0042 

674+ ' 086 
u ' D '^-0.093 

129+ - 019 
n n9q7+ - 0018 

u.uzoi _ 0013 



1.01 
0.602 



a +0.066 
-0.059 
-,+0.026 
-0.037 
,+0.049 
-0.046 





No constraint 



1.064 

0.959; 

0.592; 
< 0.84 (95%) 
No constraint 



No constraint 
0.962^ 



(-0.098 
-0.041 



1.03 



0.602 



+0.045 



0.637 



+0.12 
-0.05 

1 

+0.045 

0.046 



1.02 



0.633 



+0.16 
0.06 

1 

hO.044 
-0.041 



0.989 





+0.061 
0.031 



0.652 



+0.049 



-0.046 















No constraint No constraint No constraint 



-0.075 



-0.047 
-0.055 







-0.72 




+0.34 
-0.27 



-1 





No constraint 

-1 




fitot 

h 

T 

A s 
At 



i [Gyr] 

OS 

H 2 

H 3 

A, 

M„[CV\ 



+0.094 
0.144 
0.45 
-0.33 
0.022 
0.021 
+0.27 
0.12 
+0.17 
-0.17 

25.9±|;1 
l-14l°o: 4 3 ? 



1.095 
0.57 

0.128 
0.48 
0.33 



1.086 
0.55 

0.132 
0.50 
0.22 
20.1 
0.97 



0.14; 



-0.13 



+0.057 
0.128 
+0.47 
-0.29 
0.021 
0.028 
0.16 
0.13 
+0.34 
-0.13 
+9.2 
8.3 
+0.73 
0.23 









+0.082 
-0.058 
+0.020 
0.018 
n o 4 +0.12 

+0.13 
0.09 

17.1±|-| 
0.87«; 28 
^n +() ' 22 

u - ciu -0.17 



0.177 
0.123 



0.19] 



0.313 
0.144 
0.674 
0.15 
15.5 
0.81 





0.097 
-0.087 
0.018 
0.016 
0.087 
0.049 
+0.18 
-0.07 
+8.6 
5.6 
+0.35 
0.13 







+0.17 
-0.14 
+0.020 
0.019 
+0.14 
0.10 
+0.18 
-0.10 

94+ 40 
U - M -0.18 





0.36 
0.143 
0.63 
0.19 



0.25 
0.140 
0.74 
0.21 





0.10 
-0.10 
0.020 
0.018 
0.18 
0.07 
+0.24 
-0.11 







0.093 
-0.086 
0.020 
0.018 



I9.6t£; 8 



0.98 



+0.56 
0.21 





-,+0.26 
-0.12 



No constraint 
16.5j 
0.90j 
0.441^ 
0.424;! 
0.595 



No constraint No constraint No constraint No constraint No constraint 



+2.6 
-3.1 
0.13 
0.13 
0.013 
0.014 
+0.043 
0.040 
+0.056 
-0.048 





+2.3 
-1.8 
0.15 
0.13 
0.0090 
0.0083 
0.033 
0.029 
n «>q+0.055 



16.3 
0.87 
0.4581 
0.455 



+0.41 
0.47 
+0.17 
0.17 
+0.0067 
0.0081 
+0.034 
0.033 

n 5S4+ ' 050 



13.00 
0.84 
0.4541 
0.452 



1 ^ 7^+0.36 
lO. (O_ 5g 



0.32: 
0.426^ 











,+0.36 
-0.32 
r0.018 
-0.010 
441+0.039 
u -^ 1 -0.033 

602+ ' 045 

U.DUZ_ 046 

< 10.6 (95%) 



1 co+0.52 
io.oo_q gg 

q5+ 016 
u ' yo -0.14 
4541+0-0084 
U.^041_q 0Q85 

4?7 +0.036 
u -*' '-0.034 
n CQ1 +0.047 
°- 631 -0.045 



13 24+ ' 41 
ld - z4 -0.89 



0.99 H 



-0.19 
-0.14 
4543+0.OO83 

U.10tO_g 0Q85 
HO. 037 
-0.033 



0.474 
0.624" 



0.326 
0.153 

0.6841°;°!° 
< 0.35 (95%) 
< 25 (95%) 
0.80 


No constraint 
13.41±g;i? 
0-94±g;il 

0.0085 
0.0086 
0.032 
-0.030 
+0.048 
0.046 



0.4541 
0.475 
0.652 















X 2 /dof 



1426.1/1339 1428.4/1341 1430.9/1341 1431.8/1341 1431.8/1341 1431.5/1342 972.4/893 



t < 0.3), which cuts off much of the banana. Indeed, Fig- 
ure 2 shows that the distribution for several parameters 
(notably the reionization rcdshift z; on ) are bimodal, so 
this prior eliminates the rightmost of the two bumps. In 
the present paper, we wish to keep assumptions to a min- 
imum and therefore break the degeneracy using the SDSS 
measurements instead. Figure 5 illustrates the physical 
reason that this works so well: SDSS accurately measures 
the -P(fc) "shape parameter" T = hVt m = 0.21 ± 0.03 at 
2cr [20] , which crudely speaking determines the horizontal 
position of P{k) and this allowed region in the (f2 m , h)- 
plane intersects the CMB banana at an angle. Once E- 
polarization results from WMAP become available, they 
should provide another powerful way of breaking this de- 
generacy from WMAP alone, by directly constraining r 
— from our WMAP+SDSS analysis, we make the predic- 
tion t < 0.29 at 95% confidence for what this measure- 
ment should find. (Unless otherwise specified, we quote 
la limits in text and tables, whereas the 2-dimensional 
figures show 2a limits.) 



Figure 5 shows that the banana is well fit by h = 
0.7(rW0.3)-°- 35 , so even from WMAP+SDSS alone, we 
obtain the useful precision constraint ft(O m /0.3) ' 35 — 
0.697+°;°^ (68% ). 



B. Consistency with other measurements 

Figure 3 shows that the WMAP+SDSS allowed value 
of the baryon density uob = 0.023 ± 0.001 agrees well with 
the latest measurements u>b — 0.022 ± 0.002 from Big 
Bang Nucleosynthesis [41-43]. It is noteworthy that the 
WMAP+SDSS preferred value is higher than the BBN 
preferred value Ub = 0.019 + 0.001 of a few years ago [44], 
so the excellent agreement hinges on improved reaction 
rates in the theoretical BBN predictions [42] and a slight 
decrease in observed deuterium abundance. This is not to 
be confused with the more dramatic drop in inferred deu- 
terium abundance in preceding years as data improved, 
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Table 3: la constraints on cosmological parameters combining CMB and SDSS information. The columns compare different theoretical 
priors indicated by italics. The second last column drops the polarized WMAP information and the last column drops all WMAP 
information, replacing it by pre- WMAP CMB experiments. The 6par+«j column includes SN la information. 







; SDSS + WMAP temperature and polarization information 


No pol. 


No WMAP 




6par+f^fc + r + a 


6par+r2fc 


6par+y 


6Dar-4- fi/ 


QpcLY-\-W 


6par 


6par 


6par 




0.53+J-?? 


omt f 2 


0.776lHf 6 
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0.776I 121 


0.801°!° 


0-780l°?i9 


> 0.63 (95%) 


> 0.71 (95%) 
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600+°-°." 


5982+?-???! 

u.<jaoi_ .0032 


5948+? 1 ???? 

u.oy±o_ 0030 


5954+°-°°?! 

U ' J3J -0 .0038 


5965+° °??i 

u.o»uo_ .0030 


0.59681°-^ 


5977+0.0048 

u. t i _q Q045 




0.660+?-??? 

—0.097 


0.653+?-??? 

—0.084 


0.727+?-??i 

— 0.042 


0.620+°-°^ 

—0.087 


0.706+°°?? 

— 0.033 


0.699+?-??? 

— 0.045 


0.684+°-°!i 

—0.046 


0.691+?°?? 

—0.053 


h 2 n d 


0.103+S-^S 

— 0.022 


103+??" 

u - luJ -0 .024 


1195+?-???? 

u.x±»vj_o 0082 


135+°-°" 

o-io.J_o.012 


124+°-° 12 
"■ J - 4 — 0.011 


1222+"' 00 ?? 
u.±z,___o 0082 


0.1254+°-°^? 

— 0.0083 


0.1252+?-??? 8 

— 0.0076 


h 2 Q b 


0.0238+????? 

—0.0026 


0232+°'°°" 

u.u^o^_q 0017 


0242+?-??! 7 


0234+° °°" 

u.u^,ot_o.o011 


0232+°-°°" 
u.u_uz._o.ooiO 


0232+° °°;? 

u.uz.u__o ooiO 


0.0231+°°° n i 

—0.0009 


0.0229+?-??!? 

—0.0015 


u 
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n s 


0-97+gio 


0.98+°" 

— U.04 


1-012 + ^36 
— O.Ooo 


0.972+°°27 

— O.Uz i 


0.976+?-°?° 

— 0.0Z4 


0.977+°°?, 

— U.02b 


0.9731°;°?° 


i.oi5l?;g 3 3 6 3 


Jit + 1 


0.9852+°°°?j 


i 


0.976+°-°" 


i 


i 


l 


1 


l 
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0-584l°;°33 


0.584+?-??? 

— 0.028 


0.635+°-°?? 

— 0.021 


0.645+?-??? 

— 0.026 


0.637+?'?? 7 

— 0.027 


0.633+°°?* 

— 0. 022 


0.637l°;° 2 2 5 3 


0.588l?;° 25 5 


r 


< 0.50 (95%) 





< 0.47 (95%) 

















b 




1.03+°"" 

— 0. Id 


0.963+°-°8i 
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0.40i°' ° 


0.406+S-S2? 

— 0.091 


0.273+°-°*? 

— 0.U41 


0.380+°-°?! 

— 0.074 


0.294+°°?? 

— 0.032 


0.301+°-°!? 

— 0.042 


0.3161°-°!? 


0.309l°°f 9 




0.126±g;gig 


0.126+2-SiS 

— U.U19 


0.1438+°-°oIn 

— U.UUoU 


0.158+°-°" 

— O.Olz 


0.147+°-°; 2 

— 0.011 


.1454+°-°^ 2 


0.1486+°-°^ 

— U.UUS4 


0.1481+°-°°?! 


h 


O.SS^Q Qg 


0.550+2-222 

— O.Ooo 


0.725+°-°36 

— O.Ooo 


0.645+°-°4o 

— 0.040 


0.708+°°?? 

— O.OdO 


0.695+°°?? 

— O.Uol 


0.685l°;°? 3 8 


0.693l°;° 38 


T 
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— U. 1U 
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13.6l|| 




< 20 (95%) 


< 18 (95%) 
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At 
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t [Gyr] 
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1 ^ Q9+0.27 
j-O.o^ o 33 


13.65l°;i 


13 47+0-26 


jo c 4 +0.23 
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io kc+0.21 
lo.OO_0.23 


,o 51 +0.32 
io.oi_o 31 


C8 


o.9it :io 


0.861°;!? 
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U.OZO _ 077 


0.9281°;^ 
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fl 070+O.O88 
U.879_ 062 


fl S49+0.069 
U - S42 -0.053 


H 2 


n 441 +0-013 


D 4 t;77+0.0086 
'-0.0082 


n 4KQK+0.0081 

U.4OOO_ 0084 


4591+0.0091 


0.4545l°;° o 8 9 o 


455n+°-°° 83 


4549+0-0082 
U.4Oty_ 0083 


n 475+0-018 
U.4IO_ Q 020 


H 3 


422+ ' 027 
"■^-0.(131 


n 444+0-026 


0.468l°;g! 9 7 


0.472l°;° 22 


461+ - 018 

U.4Dl_o Q17 


459+0-018 


460+ - 017 
U.1DU_ 015 


4 8 5+ - 020 


A* 


fl 587+0-049 
U.b8C_ 041 


582+ ' 041 


0.632l°;° 2 2 ? 


0.648l°;g 2 2 8 5 


639+ - 027 
U.DOy_ 28 


n fi o 5 +0.024 
U.DoO_ 22 


639+0-024 
u.ooy_ 022 


586+ - 024 
U.OSO_ 025 


M v [eV] 











< 1.74 (95%) 














X 2 /doi 


1444.4/1357 


1445.4/1359 


1446.9/1359 


1447.3/1359 


1622.0/1531 


1447.2/1360 


987.8/911 


134.6/163 



which raised the uj\, prediction from u>b = 0.0125±0. 00125 
[45, 46]. 

The existence of dark matter could be inferred from 
CMB alone only as recently as 2001 [22], cf. [47], yet 
Figure 4 shows that WMAP alone requires dark matter 
at very high significance, refuting the suggestion of [48] 
that an alternative theory of gravity with no dark matter 
can explain CMB observations. 

Table 3 shows that once WMAP and SDSS arc com- 
bined, the constraints on three of the six vanilla parame- 
ters (wf,, LUd and n s ) are quite robust to the choice of the- 
oretical priors on the other parameters. This is because 
the CMB information that constrains them is mostly the 
relative heights of the first three acoustic peaks, which 
are left unaffected by all the other parameters except a. 
The four parameters (£lk,r,w, f v ) that are fixed by pri- 
ors in many published analyses cause only a horizontal 
shift of the peaks (Qk and w) and modified CMB power 
on larger angular scales (late ISW effect from fifc and w, 
tensor power from r). 



Figure 5 illustrates that two of the most basic cosmo- 
logical parameters, Q m and h, are not well constrained 
by WMAP alone even for vanilla models, uncertain by 
factors of about two and five, respectively (at 95% confi- 
dence). After including the SDSS information, however, 
the constraints are seen to shrink dramatically, giving 
Hubble parameter constraints h w 0.701q'q3 that are even 
tighter than (and in good agreement with) those from the 
HST project, h = 0.72 ± 0.07 [49], which is of course a 
completely independent measurement based on entirely 
different physics. (But see the next section for the cru- 
cial caveats.) Our results also agree well with those from 
the WMAP team, who obtained h w 0.73 ± 0.03 [6] by 
combining WMAP with the 2dFGRS. Indeed, our value 
for h is about 1 a lower. This is because the SDSS 
power spectrum has a slightly bluer slope than that of 
2dFGRS, favoring slightly higher f2 m -values (we obtain 
O m = 0.30 ± 0.04 as compared to the WMAP + 2dFGRS 
value Vl m — 0.26 ± 0.05). As discussed in more detail 
in Section VIII, this slight difference may be linked to 
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Table 4: la constraints on cosmological parameters as progressively more information/assumptions are added. First column uses WMAP 
data alone and treats the 9 parameters (r, Cl^, f^A> ^di w 6> A s , n s , a, r) as unknown, so the only assumptions are f v = 0, w = — 1. Moving 
to the right in the table, we add the assumptions r = a = 0, then add SDSS information, then add SN la information, then add the 
assumption that r < 0.3. The next two columns are for 6-parameter vanilla models (Q^ = r = a = 0), first using WMAP+SDSS data 
alone, then adding small-scale non-WMAP CMB data. The last two columns use WMAP+SDSS alone for 5-parameter models assuming 
n s = 1 ("vanilla lite") and n s = 0.96, r = 0.15 (V oc <f> 2 stochastic eternal inflation), respectively. 
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differences in modeling of non-linear redshift space dis- 
tortions and bias. For a thorough and up-to-date review 
of recent h— and f2 m -determinations, see [6]. 

Whereas the constraints of tub, 0Jd and n s are rather 
robust, we will see in the following section that our con- 
straints on h and £l m hinge crucially on the assump- 
tion that space is perfectly flat, and become substantially 
weaker when dropping that assumption. 

The last columns of Table 3 demonstrate excellent con- 
sistency with pre- WMAP CMB data (Appendix A. 3), 
which involves not only independent experiments but also 
partly independent physics, with much of the information 
coming from small angular scales I > 600 where WMAP 
is insensitive. In other words, our basic results and er- 
ror bars still stand even if we discard either WMAP or 
pre- WMAP data. Combining WMAP and smaller-scale 
CMB data (Table 4, 3rd last column) again reflects this 



consistency, tightening the error bars around essentially 
the same central values. 

Figure 6 compares various constraints on the linear 
clustering amplitude ag. Constraints from both galaxy 
clusters [50-52] (black) and weak gravitational lens- 
ing [53-55] (green/grey) are shown as shaded bands in 
the (O m , <78)-planc for the recent measurements listed 
in Table 5 and are seen to all be consistent with the 
WMAP+SDSS allowed region. However, we see that 
there is no part of the allowed region that simultane- 
ously matches all the cluster constraints, indicating that 
cluster-related systematic uncertainties such as the mass- 
temperature relation may still not have been fully prop- 
agated into the quoted cluster error bars. 

Comparing Figure 6 with Figure 2 from [68] demon- 
strates excellent consistency with an analysis combining 
the weak lensing data of [54] (Table 5) with WMAP, 
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Scalar spectral index n„ Dark matter density £J d 



FIG. 3: 95% constraints in the (n s ,i^b) plane. The shaded dark 
red/grey region is ruled out by WMAP alone for 6-parameter 
"vanilla" models, leaving the long degeneracy banana discussed 
in the text. The shaded light red/grey region is ruled out when 
adding SDSS information. The hatched band is required by Big 
Bang Nucleosynthesis (BBN). From right to left, the three vertical 
bands correspond to a scale-invariant Harrison-Zcl'dovich spectrum 
and to the common inflationary predictions n a = 1 — 2/N 0.96 
and n s = 1 — 3/JV ~ 0.94 (Table 6), assuming that the number 
of e-foldings between horizon exit of the observed fluctuations and 
the end of inflation is 50 < N < 60. 



small-scale CMB data and an w^-prior from Big Bang Nu- 
cleosynthesis. Figure 6 also shows good consistency with 
f2 m ~estimates from cluster baryon fractions [8], which in 
turn are larger than estimates based on mass-to-light ra- 
tio techniques reported in [8] (see [69] for a discussion of 
this). 

The constraints on the bias parameter b in Tables 3 
and 4 refer to the clustering amplitude of SDSS L* 
galaxies at the effective redshift of the survey relative 
to the clustering amplitude of dark matter at z = 0. 
If we take z ~ 0.15 as the effective redshift based on 
Figure 31 in [20], then the "vanilla lite" model (sec- 
ond last column of Table 4) gives dark matter fluctua- 
tions 0.925 times their present value and hence a phys- 
ical bias factor K = &/0-925 = 0.918/0.925 « 0.99, in 
good agreement with the completely independent mea- 
surement &* = 1.04 ± 0.11 [21] based on the bispectrum 
of L* 2dFGRS galaxies. A thorough discussion of such 
bias cross-checks is given by [70]. 



FIG. 4: 95% constraints in the (uj, uj) plane. Shaded dark 
red/grey region is ruled out by WMAP alone for 6-parameter 
"vanilla" models. The shaded light red/grey region is ruled out 
when adding SDSS information. The hatched band is required by 
Big Bang Nucleosynthesis (BBN). 



Table 5: Recent constraints in the (f2 m , <r 8 )-planc. 



Analysis 



Measurement 



Clusters: 

Voevodkin & Vikhlinin '03 [56] ag 

Bahcall & Bode '03, z < 0.2 [50] ag 

Bahcall & Bode '03, z > 0.5 [50] ag 

Pierpaoli et al. '02 [57] ag 

Allen et al. '03 [52] ag 

Schueckeret al. '02 [58] ag 

Viana et al. '02 [59] <r 8 

Seljak '02 [60] ag 

Reiprich & Bohringer '02 [61] ag 

Borgani et al. '01 [62] ag 

Pierpaoli et al. '01 [51] erg 



= 0.60 + 0.28f^ 5 ±0.04 
(fW0.3)°- 60 = 0.68 ±0.06 

(fW0.3) - 14 =0.92 ±0.09 

— n 77+0-05 

— u -' '-0.04 

(n m /0.3)° 



0.69 ±0.04 



: 0.711 
: 0.78 



(O m /0.3) 



+0.039 
0.031 
+0.15 



(for n m = 0.35) 
= 0.77 ±0.07 



: 0.96 



+0.15 



0.12 

= 0.66 ± -0.06 
(fW0.3)°- e0 = I- 11 - 



+0.070 
0.076 



Weak lensing: 

Heymans et al. '03 
Jarvis et al. '02 
Brown et al. '02 
Hoekstra et al. '02 
Refregieret al. '02 
Bacon et al. '02 
Van Waerbeke et al. '02 

Hamana et al. '02 



[63] CT8 (fW0.3) 

[64] CT8 (fW0.3) 

[53] <r 8 (fW0.3) 

[54] CT8 (n m /0.3) 

[65] <r 8 (fW0.3) 

[55] ag 

[66] <r 8 



0.6 



0.67 ± 0.10 



: 0.71 



+0.06 
0.08 



0.74 ±0.09 



: 0.86 
: 0.94 



+0.05 
0.07 
+0.24 



(n m /0.3) - 68 = 0.97 ± 0.13 

(n m /o.3)(°- 2 ^°- 18 ^- OA9 

0.94 4 



,+0.14 
-0.12 



[67] <r 8 (n m /0.3)- 



(0-78lg;%) 
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FIG. 5: 95% constraints in the (f2 m ,h) plane. Shaded dark 
red/grey region is ruled out by WMAP alone for 6-parameter 
"vanilla" models, leaving the long degeneracy banana discussed 
in the text. The shaded light red/grey region is ruled out when 
adding SDSS information, which can be understood as SDSS accu- 
rately measuring the P(k) "shape parameter" hf2 m = 0.21 ± 0.03 
at 2cr (sloping hatched band). The horizontal hatched band is re- 
quired by the HST key project [49]. The dotted line shows the fit 
h = 0.7(f2 m /0.3) — 0,35 , explaining the origin of the accurate con- 
straint /i(n m /0.3) - 35 = 0.70 ± 0.01 (la). 



IV. CURVED MODELS 

Let us now spice up the vanilla model space by adding 
spatial curvature flk as a free parameter, both to con- 
strain the curvature and to quantify how other con- 
straints get weakened when dropping the flatness as- 
sumption. 

Figures 7 and 8 show that there is a strong degeneracy 
between the curvature of the universe Qk = 1 — fitot and 
both the Hubble parameter h and the age of the universe 
to, when constrained by WMAP alone (even with only 
the seven parameters we are now considering allowed to 
change); without further information or priors, one can- 
not simultaneously demonstrate spatial flatness and mea- 
sure h or to. We see that although WMAP alone abhors 
open models, requiring f2 to t = + = 1 — ^fe ~ 0-9 
(95%), closed models with f2 tot as large as 1.4 are still 
marginally allowed provided that the Hubble parameter 
h ~ 0.3 and the age of the Universe to ~ 20 Gyr. Al- 
though most inflation models do predict space to be flat 
and closed inflation models require particularly ugly fine- 
tuning [71], a number of recent papers on other subjects 
have considered nearly-flat models either to explain the 




FIG. 6: 95% constraints in the (O m ,crg) plane. Shaded dark 
red/grey region is ruled out by WMAP alone for 6-parameter 
"vanilla" models. The shaded light red/grey region is ruled out 
when adding SDSS information. The 95% confidence regions are 
hatched for various recent cluster (black) and lensing (green/grey) 
analyses as discussed in the text. The vertical lines indicate the 
constraints described in [8] from mass-to-light ratios in galaxies 
and clusters (0.06 < Q m < 0.22) and from cluster baryon fractions 
(0.22 < Q m < 0.37). 



low CMB quadrupole [72] or for anthropic reasons [73- 
75] , so it is clearly interesting and worthwhile to test the 
flatness assumption observationally. In the same spirit, 
measuring the Hubble parameter h independently of the- 
oretical assumptions about curvature and measurements 
of galaxy distances at low redshift provides a powerful 
consistency check on our whole framework. 

Including SDSS information is seen to reduce the cur- 
vature uncertainty by about a factor of three. We also 
show the effect of adding the above-mentioned prior 
t < 0.3 and SN la information from the 172 SN la com- 
piled by [35], which is seen to further tighten the cur- 
vature constraints to n to t = 1-01 ± 0.02 (lcr), providing 
a striking vindication of the standard inflationary pre- 
diction f^tot = 1. Yet even with all these constraints, 
a strong degeneracy is seen to persist between curvature 
and h, and curvature and to, so that the HST key project 
[49] remains the most accurate measurement of ft.. If we 
add the additional assumption that space is exactly flat, 
then uncertainties shrink by factors around 3 and 4 for 
h and to, respectively, still in beautiful agreement with 
other measurements. The age limit to > 12 Gyr shown in 
Figure 8 is the 95% lower limit from white dwarf ages by 
[76]; for a thorough reviews of recent age determinations, 
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FIG. 7: 95% constraints in the (fitot,/*) plane. Shaded dark 
red/grey region is ruled out by WMAP alone for 7-parameter 
curved models, showing that CMB fluctuations alone do not simul- 
taneously show space to be flat and measure the Hubble parameter. 
The shaded light red/grey region is ruled out when adding SDSS 
information. Continuing inwards, the next two regions are ruled 
out when adding the r < 0.3 assumption and when adding SN la 
information as well. The light hatched band is required by the HST 
key project [49]. The dotted line shows the fit h = 0.7fi^, explain- 
ing the origin of the accurate constraints hQ,^ ot = 0.703^q'q24 and 
n to t(V0-7) - 2 = 1.001+°;°°* (la). 



see [6, 77]. 

This curvature degeneracy is also seen in Figure 9, 
which illustrates that the existence of dark energy Oa > 
is only required at high significance when augment- 
ing WMAP with either galaxy clustering information or 
SN la information (as also pointed out by [6]). This 
stems from the well-known geometric degeneracy where 
ilk and f?A can be altered so as to leave the acous- 
tic peak locations unchanged, which has been exhaus- 
tively discussed in the pre- WMAP literature — see, e.g., 
[12, 13, 15, 16, 78]. 

In conclusion, we obtain sharp constraints on spatial 
curvature and interesting constraints on h, to and Q\, but 
only when combining WMAP with SDSS and/or other 
data. In other words, within the class of almost flat 
models, the WMAP-only constraints on h, to and £l\ 
are weak, and including SDSS gives a huge improvement 
in precision. 

Since the constraints on h and to are further tight- 
ened by a large factor if space is exactly flat, can one 
justify the convenient assumption f2 to t = 1? Although 
WMAP alone marginally allows f2 to t = 1.5 (Figure 7), 



FIG. 8: 95% constraints in the (fitot, to) plane. Shaded dark 
red/grey region is ruled out by WMAP alone for 7-parameter 
curved models, showing that CMB fluctuations do not simulta- 
neously show space to be flat and measure the age of the Universe. 
The shaded light red/grey region is ruled out when adding SDSS 
information. Continuing inwards, the next two regions are ruled 
out when adding the r < 0.3 assumption and when adding SN la 
information as well. Stellar age determinations (see text) rule out 
t < 12 Gyr. 



WMAP+SDSS shows that fl tot is within 15% of unity. 
It may therefore be possible to bolster the case for per- 
fect spatial flatness by demolishing competing theoretical 
explanations of the observed approximate flatness — for 
instance, it has been argued that if the near-flatness is 
due to an anthropic selection effect, then one expects de- 
partures from fitot ~ 1 of order unity [73, 75], perhaps 
larger than we now observe. This approach is particu- 
larly promising if one uses a prior on h. Imposing a hard 
limit 0.58 < h < 0.86 corresponding to the 2er range from 
the HST key project [49], we obtain O tot = 1.030±g;g|| 
from WMAP alone, fi tot = 1.023t{J;{g§ adding SDSS 
and f^tot = l-OlO^o oi7 when also adding SN la and the 
t < 0.3 prior. 1 



Within the framework of Bayesian inference, such an argument 
would run as in the following example. Let us take the current 
best measurement from above to be fitot = 1.01 ±0.02 and use it 
to compare an inflation model predicting fitot = 1 ± 10 — 5 with a 
non-inflationary FRW model predicting that a typical observer 
sees fitot = 1 ± 1 because of anthropic selection effects [73—75]. 
Convolving with the 0.02 measurement uncertainty, our two rival 
models thus predict that our observed best-fit value is drawn 
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FIG. 9: 95% constraints in the (Q m ,Cl\) plane. Shaded dark 
red/grey region is ruled out by WMAP alone for 7-parameter 
curved models, illustrating the well-known geometric degeneracy 
between models that all have the same acoustic peak locations. 
The shaded light red/grey region is ruled out when adding SDSS 
information. Continuing inwards, the next two regions are ruled 
out when adding the r < 0.3 assumption and when including SN 
la information as well. Models on the diagonal dotted line are flat, 
those below are open and those above are closed. The constraints 
in this plot agrees well with those in Figure 13 from [6] when taking 
the t prior into account. 



V. TESTING INFLATION 

A. The generic predictions 

Two generic predictions from inflation are perfect flat- 
ness (flk = 0, i.e., f^tot = 1 — f2fc = 1) and approxi- 
mate scale-invariance of the primordial power spectrum 



from distributions fitot = li0.02 and Ctot = 1±1, respectively. 
If we approximate these distributions by Gaussians /(Htot) = 
e -m tot -i)/<r] 2 /2/^2^ a . with a _ q 02 anc i a = 1, respectively, 
we find that the observed value is about 22 times more likely 
given inflation. In other words, if we view both models as equally 
likely from the outset, the standard Bayesian calculation 



Explanation 


Prior prob. 


Obs. 


likelihood 


Posterior prob. 


Inflation 


0.5 




17.6 


0.96 


Anthropic 


0.5 




0.80 


0.04 



strongly favors the inflationary model. Note that it did not have 
to come out this way: observing Qtot = 0.90 ± 0.02 would have 
given 99.99% posterior probability for the anthropic model. 



(n s ~ 1). Tables 2-4 show that despite ever-improving 
data, inflation still passes both of these tests with flying 
colors. 2 

The tables show that although all cases we have ex- 
plored are consistent with f2 t ot = n s = 1, adding priors 
and non-CMB information shrinks the error bars by fac- 
tors around 6 and 4 for f2 to t and n s , respectively. 

For the flatness test, Table 4 shows that S7 to t is within 
about 20% of unity with 68% confidence from WMAP 
alone without priors (even f2 t ot ~ 1-5 is allowed at the 
95% confidence contour). When we include SDSS, the 
68% uncertainty tightens to 10%, and the errors shrink 
impressively to the percent level with more data and pri- 
ors: n to t = l-012t° 022 usm S WMAP, SDSS, SN la and 
r < 0.3. 

For the scalar spectral index, Table 4 shows that n s ~ 1 
to within about 15% from WMAP alone without priors, 
tightening to n s = 0.977^^025 when adding SDSS and as- 
suming the vanilla scenario, so the cosmology community 
is rapidly approaching the milestone where the depar- 
tures from scale-invariance that most popular inflation 
models predict become detectable. 



B. Tensor fluctuation 

The first really interesting confrontation between the- 
ory and observation was predicted to occur in the (n s , r) 
plane (Figure 10), and the first skirmishes have already 
begun. The standard classification of slow-roll inflation 
models [80-82] characterized by a single field inflaton po- 
tential V(4>) conveniently partitions this plane into three 
parts (Figure 10) depending on the shape of V(4>): 

1. Small-field models are of the form expected from 
spontaneous symmetry breaking, where the poten- 
tial has negative curvature V{(j>)" < and the field 
<fi rolls down from near the maximum, and all pre- 
dict r < |(1 — n s ), n s < 1. 

2. Large-field models are characteristic of so-called 
chaotic initial conditions, in which <f> starts out far 
from the minimum of a potential with positive cur- 
vature {y"(4>) > 0), and all predict §(1 - n s ) < 
r < 8(1 - n s ), n s < 1. 

3. Hybrid models are characterized by a field rolling 
toward a minimum with V =/= 0. Although they 
generally involve more than one inflaton field, they 
can be treated during the inflationary epoch as 
single-field inflation with V" > and predict r > 
|(1 —n s ), also allowing n s > 1. 



2 Further successes, emphasized by the WMAP team and [79], 
are the inflationary predictions of adiabaticity and phase coher- 
ence which account for the peak/trough structure observed in 
the CMB power spectrum. 
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Scalar spectral index n s 

FIG. 10: 95% constraints in the (n s ,r) plane. Shaded dark 
red/grey region is ruled out by WMAP alone for 7-parameter mod- 
els (the vanilla models plus r). The shaded light red/grey region 
is ruled out when adding SDSS information. The two dotted lines 
delimit the three classes of inflation models known as small-field, 
large-field and hybrid models. Some single-field inflation models 
make highly specific predictions in this plane as indicated. Prom 
top to bottom, the figure shows the predictions for V(rp) oc </> 6 
(line segment; ruled out by CMB alone), V(<j>) oc <f> 4 (star; a text- 
book inflation model; on verge of exclusion) and V(<f>) oc <f> 2 (line 
segment; the eternal stochastic inflation model; still allowed), and 
V(<f>) oc 1 — (4>/<f>*) 2 (horizontal line segment with r ~ 0; still al- 
lowed). These predictions assume that the number of e- foldings 
between horizon exit of the observed fluctuations and the end of 
inflation is 64 for the <f> 4 model and between 50 and 60 for the 
others as per [86]. 



These model classes are summarized in Table 6 together 
with a sample of special cases. For details and deriva- 
tions of the tabulated constraints, see [5, 80-85]. For 
comparison with other papers, remember that we use the 
same normalization convention for r as CMBfast and the 
WMAP team, where r = — 8nt for slow-roll models. The 
limiting case between small-field and large-field models 
is the linear potential V((f>) oc (f>, and the limiting case 
between large-field and hybrid models is the exponential 
potential V(<j>) oc e^*. The WMAP team [5] further 
refine this classification by splitting the hybrid class into 
two: models with n s < 1 and models with n s > 1. 

Many inflationary theorists had hoped that early data 
would help distinguish between these classes of models, 
but Figure 10 shows that all three classes are still allowed. 

What about constraints on specific inflation models as 
opposed to entire classes? Here the situation is more 
interesting. Some models, such as hybrid ones, allow 



two-dimensional regions in this plane. Table 6 shows 
that many other models predict a one-dimensional line or 
curve in this plane. Finally, a handful of models are ex- 
tremely testable, making firm predictions for both n s and 
r in terms of N, the number of e-foldings between horizon 
exit of the observed fluctuations and the end of inflation. 
Recent work [86, 87] has shown that 50 < N < 60 is re- 
quired for typical inflation models. The quartic model 
V ^ cf) 4 is an anomaly, requiring N w 64 with very 
small uncertainty. Figure 10 shows that power law mod- 
els V oc <j) p are ruled out by CMB alone for p = 6 and 
above. Figure 10 indicates that the textbook V oc tp 4 
model (indicated by a star in the figure) is marginally 
allowed. [5] found it marginally ruled out, but this as- 
sumed N — 50 — the subsequent result N « 64 [86] 
pushes the model down to the right and make it less dis- 
favored. V oc (j> 2 has been argued to be the most natural 
power-law model, since the Taylor expansion of a generic 
function near its minimum has this shape and since there 
is no need to explain why quantum corrections have not 
generated a quadratic term. This potential is used in the 
stochastic eternal inflation model [88], and is seen to be 
firmly in the allowed region, as are the small-field "tomb- 
stone model" from Table 6 and the GUT-scale model of 
[89] (predicting n s = 1 - 1/N « 0.98, r w 10" 8 ). 

In conclusion, Figure 10 shows that observations are 
now beginning to place interesting constraints on infla- 
tion models in the (n s , r)-plane. As these constraints 
tighten in coming years, they will allow us to distin- 
guish between many of the prime contenders. For in- 
stance, the stochastic eternal inflation model predicting 
(n s ,r) w (0.96,0.15) will become distinguishable from 
models with negligible tensors, and in the latter cate- 
gory, small-field models with, say, n s < 0.95, will become 
distinguishable from the scale- invariant case n s = 1. 



C. A running spectral index? 

Typical slow-roll models predict not only negligible 
spatial curvature, but also that the running of the spec- 
tral index a is unobservably small. We therefore assumed 
f2fe = a = when testing such models above. 

Let us now turn to the issue of searching for departures 
from a power law primordial power spectrum. This is- 
sue has generated recent interest after the WMAP team 
claim that a < was favored over a = 0, at least at 
modest statistical significance, with the preferred value 
being a ~ —0.07 [5, 6]. 

Slow-roll models typically predict \a\ of order A -2 ; 
for these models, |a| is rarely above 10 -3 , much smaller 
than the WMAP-team preferred value. Those inflation 
models that do predict such a strong second derivative of 
the primordial power spectrum (in log-log space) tend to 
produce substantial third and higher derivatives as well, 
so that a parabolic curve parametrized by A s , n s and a 
is a poor approximation of the model (e.g., [90]). Lack- 
ing strong theoretical guidance one way or another, we 
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Table 6: Sample inflation model predictions. N is the number of e-folds between horizon exit of the observed fluctuations and the end of 

inflation. 
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therefore drop our priors on flk and r when constraining 
a. 

Tables 2 and 3 show that our best-fit a-values agree 
with those of [5], but are consistent with a — 0, since 
the 95% error bars are of order 0.1. They show that x 2 
drops by only 5 relative to vanilla models, which is not 
statistically significant because a drop of 3 is expected 
from freeing the three parameters f^, r and a. More- 
over, we see that our WMAP-only constraint is similar 
to our WMAP+SDSS constraint, showing that any hint 
of running comes from the CMB alone, most likely from 
the low quadrupole power [6]; see also [91, 92]. This is 
at least qualitatively consistent with the WMAP team 
analysis [6]; apart from the low quadrupole, most of the 
evidence that a ^ comes from CMB fluctuation data on 
small scales (e.g., the CBI data [93]) and measurements 
of the small-scale fluctuations from the Lya forest; in- 
deed, including the 2dFGRS data slightly weakens the 
case for running. For the Lya forest case, the key issue is 
the extent to which the measurement uncertainties have 
been adequately modeled [94], and this should be clar- 
ified by the forthcoming Lya forest measurements from 
the SDSS. 
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VI. NEUTRINO MASS 

It has long been known [95] that galaxy surveys are 
sensitive probes of neutrino mass, since they can detect 
the suppression of small-scale power caused by neutri- 
nos streaming out of dark matter overdensities. For de- 
tailed discussion of post- WMAP astrophysical neutrino 
constraints, see [6, 96-99], and for an up-to-date review 
of the theoretical and experimental situation, see [100]. 

Our neutrino mass constraints are shown in the M v - 
panel of Figure 2, where we allow our standard 6 "vanilla" 
parameters and /„ to be free. The most favored value is 
M v = 0, and obtain a 95% upper limit M v < 1.7eV. Fig- 
ure 11 shows that WMAP alone tells us nothing whatso- 



FIG. 11: 95% constraints in the (uj^.f^) plane. Shaded dark 
red/grey region is ruled out by WMAP alone when neutrino mass 
is added to the 6 "vanilla" models. The shaded light red/grey 
region is ruled out when adding SDSS information. The five curves 
correspond to M„, the sum of the neutrino masses, equaling 1, 2, 
3, 4 and 5 eV, respectively — barring sterile neutrinos, no neutrino 
can have a mass exceeding ~ My/3. 



ever about neutrino masses and is consistent with neu- 
trinos making up 100% of the dark matter. Rather, the 
power of WMAP is that it constrains other parameters 
so strongly that it enables large-scale structure data to 
measure the small-scale P(fc)-suppression that massive 
neutrinos cause. 

The sum of the three neutrino masses (assuming stan- 
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dard freezeout) is [38] M v w (94.4 eV)u> d f v . The neu- 
trino energy density must be very close to the standard 
freezeout density [101-103], given the large mixing angle 
solution to the solar neutrino problem and near maxi- 
mal mixing from atmospheric results — see [104, 106] for 
up-to-date reviews. Any substantial asymmetries in neu- 
trino density from the standard value would be "equili- 
brated" and produce a primordial 4 He abundance incon- b 
sistent with that observed. 

Our upper limit is complemented by the lower limit ~ 
from neutrino oscillation experiments. Atmospheric neu- 6 
trino oscillations show that there is at least one neu- cao 
trino (presumably mostly a linear combination of and -| 
v T ) whose mass exceeds a lower limit around 0.05 eV *» 
[100, 104]. Thus the atmospheric neutrino data corre- ~ 
sponds to a lower limit w v > 0.0005, or /„ > 0.004. The 
solar neutrino oscillations occur at a still smaller mass 
scale, perhaps around 0.008 eV [100, 105, 106]. These 
mass-splittings are much smaller than 1.7 eV, suggest- 
ing that all three mass eigenstates would need to be al- 
most degenerate for neutrinos to weigh in near our upper 
limit. Since sterile neutrinos are disfavored from being 
thermalized in the early universe [107, 108], it can be as- 
sumed that only three neutrino flavors are present in the 
neutrino background; this means that none of the three 
neutrinos can weigh more than about 1.7/3 = 0.6 eV. 
The mass of the heaviest neutrino is thus in the range 
0.04 - 0.6 eV. 

A caveat about non-standard neutrinos is in order. 
To first order, our cosmological constraint probes only 
the mass density of neutrinos, p„, which determines the 
small-scale power suppression factor, and the velocity dis- 
persion, which determines the scale below which the sup- 
pression occurs. For the low mass range we have dis- 
cussed, the neutrino velocities are high and the suppres- 
sion occurs on all scales where SDSS is highly sensitive. 
We thus measure only the neutrino mass density, and 
our conversion of this into a limit on the mass sum as- 
sumes that the neutrino number density is known and 
given by the standard model freezeout calculation, 112 
cm~ 3 . In more general scenarios with sterile or other- 
wise non-standard neutrinos where the freezeout abun- 
dance is different, the conclusion to take away is an 
upper limit on the total light neutrino mass density of 
p v < 4.8 x 10~ 28 kg/m 3 (95%). To test arbitrary non- 
standard models, a future challenge will be to indepen- 
dently measure both the mass density and the velocity 
dispersion, and check whether they are both consistent 
with the same value of M v . 

The WMAP team obtains the constraint M u < 0.7 eV 
[6] by combining WMAP with the 2dFGRS. This limit is 
a factor of three lower than ours because of their stronger 
priors, most importantly that on galaxy bias b deter- 
mined using a bispectrum analysis of the 2dF galaxy 
clustering data[21]. This bias was measured on scales 
k ~ 0.2 — 0.4/i/Mpc and assumed to be the same on the 
scales k < 0.2/i/Mpc that were used in the analysis. In 
this paper, we prefer not to include such a prior. Since 
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FIG. 12: Constraints in the (f v ,crs) plane. Shaded dark red/grey 
region is ruled out by WMAP alone (95%) when neutrino mass is 
added to the 6 "vanilla" models. The shaded light red/grey region 
is ruled out when adding SDSS information. The recent claim that 
f v >0 [109] hinges on assuming that galaxy clusters require low erg- 
values (shaded horizontal band) and dissolves when using what we 
argue are more reasonable uncertainties in the cluster constraints. 



the bias is marginalized over, our SDSS neutrino con- 
straints come not from the amplitude of the power spec- 
trum, only from its shape. This of course allows us to 
constrain b from WMAP+SDSS directly; we find values 
consistent with unity (for L* galaxies) in almost all cases 
(Tables 3 and 4) . A powerful consistency test is that our 
corresponding value = 0.54±°;^ from WMAP+SDSS 
agrees well with the value ~ 0.5 measured from redshift 
space distortions in [20]. 

Seemingly minor assumptions can make a crucial dif- 
ference for neutrino conclusions, as discussed in detail 
in [6, 96, 97]. A case in point is a recent claim that 
nonzero neutrino mass has been detected by combining 
WMAP, 2dFGRS and galaxy cluster data [109]. Figure 
2 in that paper (middle left panel) shows that nonzero 
neutrino mass is strongly disfavored only when including 
data on X-ray cluster abundance, which is seen (lower 
middle panel) to prefer a low normalization of order 
cr 8 « 0.70 ± 0.05 (68%). Figure 12 provides intuition 
for the physical origin on the claimed neutrino mass de- 
tection. Since WMAP fixes the normalization at early 
times before neutrinos have had their suppressing effect, 
we see that the WMAP-allowed o^-value drops as the 
neutrino fraction f v increases. A very low ag-value there- 
fore requires a nonzero neutrino fraction. The particular 
cluster analysis used by [109] happens to give one of the 
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lowest cg-values in the recent literature. Table 5 and Fig- 
ure 6 show a range of values larger than the individual 
quoted errors, implying the existence of significant sys- 
tematic effects. If we expand the error bars on the cluster 
constraints to as — 0.8 ± 0.2, to reflect the spread in the 
recent literature, we find that the evidence for a cosmo- 
logical neutrino mass detection goes away. The sensitiv- 
ity of neutrino conclusions to cluster a8 normalization 
uncertainties was also discussed in [109]. 



VII. DARK ENERGY EQUATION OF STATE 

Although we now know its present density fairly accu- 
rately, we know precious little else about the dark energy, 
and post-WMAP research is focusing on understanding 
its nature [110-118]. Above we have assumed that the 
dark energy behaves as a cosmological constant with its 
density independent of time, i.e., that its equation of 
state w = — 1. Figure 2 and Figure 13 show our con- 
straints on w, assuming that the dark energy is homoge- 
neous, i.e., does not cluster 3 . Although our analysis adds 
improved galaxy and SN la data to that of the WMAP 
team [6] and uses different assumptions, Figure 13 agrees 
well with Figure 1 1 from [6] and our conclusions are qual- 
itatively the same: adding roasa free parameter does not 
help improve \ 2 f° r the best fit, and all data are consis- 
tent with the vanilla case w — —1, with uncertainties in w 
at the 20% level. [110, 117] obtained similar constraints 
with different data and an /i-prior. 

Tables 2 and 3 show the effect of dropping the w = — 1 
assumption on other parameter constraints. These ef- 
fects are seen to be similar to those of dropping the flat- 
ness assumption, but weaker, which is easy to understand 
physically. As long as there are no spatial fluctuations 
in the dark energy (as we have assumed), changing w 
has only two effects on the CMB: it shifts the acoustic 
peaks sideways by altering the angle-distance relation, 
and it modifies the late Integrated Sachs- Wolfe (ISW) ef- 
fect. Its only effect effect on the matter power spectrum 
is to change its amplitude via the linear growth factor. 
The exact same things can be said about the parameters 
and f2fc, so the angle-diameter degeneracy becomes 
a two-dimensional surface in the three-dimensional space 
(fife, f^A, w), broken only by the late ISW effect. Since 
the peak-shifting is weaker for w than for Q k (for changes 



Dark energy clustering can create important modifications of the 
CMB power spectrum and can weaken the to-constraints by in- 
creasing degeneracies [117]. We have ignored the effect of dark 
energy clustering since it depends on the dark energy sound 
speed, which is in turn model-dependent and at present com- 
pletely unknown. Indeed, all evidence for dark energy so far 
traces back to the observed cosmic expansion history H(z) de- 
parting from f2j\ = Friedmann equation, and if this departure 
is caused by modified gravity rather than some sort of new sub- 
stance, then there may be no dark energy fluctuations at all. 




FIG. 13: 95% constraints in the (Q m ,w) plane. Shaded dark 
red/grey region is ruled out by WMAP alone when equation of 
state w is added to the 6 "vanilla" parameters. The shaded light 
red/grey region is ruled out when adding SDSS information, and 
the yellow/very light grey region is excluded when including SN la 
information as well. 



generating comparable late ISW modification), adding w 
to vanilla models wreaks less havoc with, say, h than does 
adding to vanilla models (Section IV). 



VIII. DISCUSSION AND CONCLUSIONS 

We have measured cosmological parameters using 
the three-dimensional power spectrum P(k) from over 
200,000 galaxies in the Sloan Digital Sky Survey (SDSS) 
in combination with WMAP and other data. Let us first 
discuss what we have and have not learned about cosmo- 
logical parameters, then summarize what we have and 
have not learned about the underlying physics. 



A. The best fit model 

All data we have considered are consistent with a 
"vanilla" flat adiabatic ACDM model with no tilt, run- 
ning tilt, tensor fluctuations, spatial curvature or mas- 
sive neutrinos. Readers wishing to choose a concordance 
model for a calculational purposes using Ockham's razor 
can adopt the best fit "vanilla lite" model 

(r, Q A , w d , u b , A s ) = (0.17, 0.72, 0.12, 0.024, 0.89) (2) 
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(Table 4, second last column). Note that this is even 
simpler than 6-parameter vanilla models, since it has 
n s = 1 and only 5 free parameters [119]. A more the- 
oretically motivated 5-parameter model is that of the ar- 
guably most testable inflation model, V oc 4> 2 stochas- 
tic eternal inflation, which predicts (n s ,r) = (0.15,0.96) 
(Figure 10) and prefers 

(t, Oy, u d , uj b , A s ) = (0.09, 0.68, 0.123, 0.023, 0.75) (3) 

(Table 4, second last column). 

Note that these numbers are in substantial agreement 
with the results of the WMAP team [6], despite a com- 
pletely independent analysis and independent redshift 
survey data; this is a powerful confirmation of their re- 
sults and the emerging standard model of cosmology. 
Equally impressive is the fact that we get similar results 
and error bars when replacing WMAP by the combined 
pre- WMAP CMB data (compare the last columns of Ta- 
ble 3). In other words, the concordance model and the 
tight constraints on its parameters are no longer depen- 
dent on any one data set — everything still stands even if 
we discard either WMAP or pre- WMAP CMB data and 
either SDSS or 2dFGRS galaxy data. No single data set 
is indispensable. 

As emphasized by the WMAP team, it is remarkable 
that such a broad range of data are describable by such a 
small number of parameters. Indeed, as is apparent from 
Tables 2-4, \ 2 does not improve significantly upon the 
addition of further parameters for any set of data. How- 
ever, the "vanilla lite" model is not a complete and self- 
consistent description of modern cosmology; for example, 
it ignores the well-motivated inflationary arguments for 
expecting n s ^ 1. 

B. Robustness to physical assumptions 

On the other hand, the same criticism can be leveled 
against 6-parameter vanilla models, since they assume 
r = even though some of the most popular inflation 
models predict a significant tensor mode contribution. 
Fortunately, Table 3 shows that augmenting vanilla mod- 
els with tensor modes has little effect on other parame- 
ters and their uncertainties, mainly just raising the best 
fit spectral index n s from 0.98 to 1.01. 

Another common assumption is that the neutrino den- 
sity f v is negligible, yet we know experimentally that 
/„ > and there is an anthropic argument for why neu- 
trinos should make a small but non-negligible contribu- 
tion [120]. The addition of neutrinos changes the slope of 
the power spectrum on small scales; in particular, when 
we allow /„ to be a free parameter, the value of as drops 
by 10% and Q TO increases by 25% (Table 2). 

We found that the assumption with the most striking 
implications is that of perfect spatial flatness, fi to t = 1 
dropping it dramatically weakens the limits on the 
Hubble parameter and the age of the Universe, allowing 
h = 0.5 and t = 18 Gyr. Fortunately, this flatness 



assumption is well-motivated by inflation theory; while 
anthropic explanations exist for the near flatness, they 
do not predict the Universe to be quite as flat as it is 
now observed to be. 

Constraints on other parameters are also somewhat 
weakened by allowing a running spectral index a ^= 
and an equation of state w ^ — 1, but we have argued 
that these results are more difficult to take seriously 
theoretically. It is certainly worthwhile testing whether 
n s depends on k and whether Q\ depends on z, but 
parametrizing such departures in terms of constants a 
and w to quantify the degeneracy with other parameters 
is unconvincing, since most inflation models predict ob- 
servably large |a| to depend strongly on k and observably 
large \w + 1| can depend strongly on z. 

It is important to parametrize and constrain possible 
departures the current cosmological framework: any test 
that could have falsified it, but did not, bolsters its cred- 
ibility. Post- WMAP work in this spirit has included con- 
straints on the dark energy sound speed [117] and time 
dependence [121, 122], the fine structure constant [123], 
the primordial helium abundance [124, 125], isocurvature 
modes [126] and features in the primordial power spec- 
trum [127, 128]. 



C. Robustness to data details 

How robust are our cosmological parameter measure- 
ments to the choice of data and to our modeling thereof? 

For the CMB, most of the statistical power comes from 
the unpolarized WMAP data, which we confirmed by 
repeating our 6-parameter analysis without polarization 
information. The main effect of adding the polarized 
WMAP data is to give a positive detection of r (Sec- 
tion VIII D 4 below). The quantity as&~ T determines the 
amplitudes of acoustic peak amplitudes, so the positive 
detection of r leads to a value of a s 15% higher than 
without the polarization data included. 

For the galaxy P(k) data, there are options both for 
what data set to use and how to model it. To get a 
feeling for the quantitative importance of choices, we re- 
peat a simple benchmark analysis for a variety of cases. 
Let us measure the matter density f2 m using galaxy data 
alone, treating A s as a second free parameter and fixing 
all others at the values fife = f v = a = 0, tub = 0.024, 
n s = 1, w = — 1, b = 1 and h = 0.72. Roughly speaking, 
we are thus fitting the measured galaxy power spectrum 
to a power spectrum curve that we can shift horizontally 
(with our "shape parameter" fl m ) and vertically (with 
A s ). We have chosen this particular example because, 
as described in Section III, it is primarily this shape pa- 
rameter measurement that breaks the WMAP banana 
degeneracy. The parameters r and r of course have no 
effect on P(k), and the remaining two are determined 
by the matter density via the identities Oa = 1 — fi m , 

LOd = h 2 fl m — LOb. 

Our results are summarized in Table 7. We stress 
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Tabic 7: Robustness to data and method details. 



Analysis 



Baseline 

U = 0.15/i/Mpc 
U = 0.1/i/Mpc 
No bias correction 
Linear P(k) 
2dFGRS 



2gl +0.033 

u ' z31 -0 .027 
q 2g7 +0.038 
u - zyl -0.032 
n 001+0.079 
u ' JJi -0.051 

25fi+ ' 027 

U.ZOD_ q 24 
q oo 4 +0.027 

n oki +0 036 
u - zol -0.027 
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that they should not be interpreted as realistic measure- 
ments of f2 m , since the other parameters have not been 
marginalized over. This is why the error bars are seen to 
be smaller even than when WMAP was included above 
(last column of Table 4). 

To avoid uncertainties associated with nonlinear red- 
shift space distortions and scale-dependent galaxy bias, 
we have used SDSS measurements of P(k) only for k < 
fcmax throughout this paper, chosing fc max = 0.2/i/Mpc as 
recommended in [20] . The WMAP team made this same 
choice fc max = 0.2/i/Mpc when analyzing the 2dFGRS 
data [7]. An option would be to tighten this cut to be 
still more cautious. Table 7 shows that cutting back to 
J; max = 0.15/i/Mpc has essentially no effect on the best-fit 
f2 m -value and increases error bars by about 20%. Cut- 
ting back all the way down to fcmax = 0.1/i/Mpc is seen 
to more than double the baseline error bars, the baseline 
measurement lying about 0.6a below the new best fit. 

As described in [20] , the SDSS measurements were cor- 
rected for luminosity-dependent bias. Table 7 shows that 
if this were not done, Q m would drop by about 0.03, or 
lo". This correction is of course not optional. However, if 
the correction itself were somehow inaccurate at say the 
10% level, one would expect a bias in f2 m around 0.003. 

Just like the WMAP team [6, 7] , we have used the non- 
linear matter power spectrum for all our analysis. Table 7 
shows that if we had used the linear spectrum instead, 
then Q m would rise by about 0.04, or 1.3a. This happens 
because the linear power spectrum is redder, with less 
small-scale power, which can be roughly offset by rais- 
ing Q m and hence shifting the curve to the right. Like 
the above-mentioned correction for luminosity-dependent 
bias, correction for nonlinearities must be included. How- 
ever, given the large uncertainties about how biasing be- 
haves in this quasilinear regime, it may well be that this 
correction is only accurate to 25%, say, in which case we 
would expect an additional uncertainty in f2 m at the 0.01 
level. 

Finally, we have repeated the analysis using an en- 
tirely different dataset, the P(fc)-measurement from the 
2dFGRS team [11]. Although the WMAP team used 
fcmax = 0.2/i/Mpc, we used the data available online with 
fcmax = 0.15/i/Mpc here as recommended by the 2dF- 
GRS team [11]. Table 7 shows that 2dFGRS measures a 
slightly redder power spectrum than SDSS, correspond- 
ing to Q, m down by 0.04, or 1.3c. 

In conclusion, we see that a number of issues related to 
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FIG. 14: Effect of increasing the amount of SDSS data used, given 
by the maximum fc-value used. Top panel shows how the relative 
errors on various parameters shrink as more data is included. For 
the neutrino density Q v = fvfldi the absolute rather than relative 
error is shown. Bottom panel shows the ratio of systematic errors 
to statistical errors (from top panel) grows as smaller scales are 
included. This is for the extreme case where nonlinear corrections 
are present but completely ignored, which we view as a worst-case 
scenario. 



data selection and modeling can have noticeable effects 
on the results. Internally to SDSS, such effects could 
easily change ft m by as much as 0.01, and the 2dFGRS 
difference is about 0.04, or one standard deviation - 
roughly what one would expect with two completely in- 
dependent data sets. 

To quantify the effect of systematic uncertainties when 
both other parameters and WMAP data are included, 
we carry out a second testing exercise. Using the Fisher- 
matrix technique of [129], we compute how our best-fit 
parameter values shift in response to a systematic bias 
in the theoretically computed power spectrum P(k). To 
be conservative, we make the rather extreme assump- 
tion that the measurements correspond to the nonlin- 
ear power spectrum but that the analysis ignores non- 
linear corrections entirely, simply fitting to the linear 
power spectrum. Although we view this as a worst-case 
scenario, it provides an instructive illustration of how 
problems related to nonlinear redshift space distortions 
and scale-dependent biasing might scale with fc max , the 
largest fc-band included. 

Our results are shown in Figure 14. The upper panel 
shows how the constraints from WMAP alone (on left 
side of figure) gradually improve as more SDSS data are 
included. The dramatic neutrino improvement seen at 
small fcmax is due to WMAP alone leaving the neutrino 
fraction unconstrained. The other parameters where 
SDSS helps the most are seen to be to m and h, which 
can be understood based on our discussion in Section III. 
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The SDSS power spectrum we have used does not probe 
to scales much smaller than k ~ 0.2, which is why little 
further improvement is seen beyond this value. 

The lower panel shows the ratio of the above- 
mentioned systematic error to the statistical error for 
each parameter. We see that the most sensitive param- 
eter is Ct> m , which justifies our singling it out for special 
scrutiny above in Table 7 (where u> m is equivalent to Q m 
since we kept h fixed). Although n s also partially mim- 
ics the nonlinear correction and perhaps scale-dependent 
bias, it is seen to be somewhat less sensitive. Our Fisher 
matrix estimate is seen to be somewhat overly pessimistic 
for Q m , predicting that neglecting nonlinearities shifts 
w m by of order 2er for fc max = 0.2/i/Mpc whereas the 
brute force analysis in Table 7 shows the shift to be 
only about half as large even when WMAP is ignored. 
The sensitivity to h is linked to the f2 m -sensitivity by 
the banana in Figure 5. The f2„-sensitivity comes from 
the small-scale neutrino P(fc)-supprcssion being similar 
to the suppression in going from nonlinear to linear P(k)- 
modeling. 

In conclusion, as long as errors in the modeling of non- 
linear redshift distortions and bias arc not larger than 
the nonlinear correction itself, we expect our uncertain- 
ties with fc max = 0.2/i/Mpc to be dominated by statisti- 
cal rather than systematic errors. The fact that cutting 
back to fc max = 0.15/i/Mpc left our results virtually un- 
changed (Table 7) supports this optimistic conclusion. 
Indeed, Figure 14 shows that with fc max = 0.15/i/Mpc, 
the reader wanting to perform a simple analysis can even 
use the linear P(k) to good approximation. 

However, both statistical errors and the systematic er- 
rors we have discussed in this section are dwarfed by the 
effects of changing theoretical priors. For instance, Ta- 
ble 2 shows that fl m increases by 0.08 when dropping 
either the assumption of negligible neutrinos or the as- 
sumption of negligible curvature. Moreover, to place this 
in perspective, all Bayesian analysis using Monte Carlo 
Markov Chains implicitly assumes a uniform prior on 
the space of the parameters where the algorithm jumps 
around, and different authors make different choices for 
these parameters, which can make a substantial differ- 
ence. 4 



4 For instance, we use A p = A a e~ 2r where the WMAP team uses 
A a [7], and both we and the WMAP team use the CMB peak 
location parameter S where many other groups use The 
difference between these implicit priors is given by the Jacobian 
of the transformation, which describes how the volume element 
changes and generically will have variations of order unity when 
a parameter varies by a factor of two. For a parameter p that is 
tightly constrained with a small relative error |Ap/p| <C 1, this 
Jacobian becomes irrelevant. For weakly constrained parameters 
like r, however, this can easily shift the best-fit value by la. 
For example, changing to a uniform prior on the rcionization 
redshift Z[ on oc r 2 / 3 as done by [68] corresponds to using a r- 
prior oc t -1 / 3 , which strongly weights the results towards low 

T. 



A final source of potential uncertainties involves bugs 
and algorithmic errors in the analysis software. To guard 
against this, we performed two completely independent 
analyses for many of the parameter spaces that we have 
tabulated, one using the Monte Carlo Markov Chain 
method described in Appendix A (coded up from scratch) 
and the second using the publicly available CosmoMC 
package [32] with appropriate modifications. We found 
excellent agreement between the two sets of results, with 
all differences much smaller than the statistical errors 
and prior-related uncertainties. 



D. What have we learned about physics? 

The fact that any simple model fits such accurate and 
diverse measurements is impressive evidence that the ba- 
sic theoretical framework of modern cosmology is correct, 
and that we need to take its implications seriously how- 
ever surprising they may be. What are these implica- 
tions? 



1 . Inflation 

The two generic predictions of perfect flatness (|f2fc| ^ 
10~ 5 ) and near scale-invariance have passed yet another 
test with flying colors. We find no evidence for running 
tilt. We also find no evidence for gravitational waves, and 
are therefore unable to measure the tensor spectral index 
and test the inflationary consistency relation r = —8n t . 
The most interesting confrontation between theory and 
observation is now occurring in the (n s ,r) plane (Fig- 
ure 10). We confirm the conclusion [5] that most popular 
models are still allowed, notably even stochastic eternal 
inflation with its prediction (n s ,r) « (0.96,0.16), but 
modest data improvements over the next few years could 
decimate the list of viable inflationary candidates and 
rival models [130]. 



2. Dark energy 

Since its existence is now supported by three indepen- 
dent lines of evidence (SN la, power spectrum analysis 
such as ours, the late ISW effect [131-136]) and its cur- 
rent density is well known (the last column off Table 2 
gives fl\ = 0.70 ± 0.04), the next challenge is clearly 
to measure whether its density changes with time. Al- 
though our analysis adds improved galaxy and SN la data 
to that of the WMAP team [6] , our conclusions are quali- 
tatively the same: all data are consistent with the density 
being time-independent as for a simple cosmological con- 
stant (w = —1), with uncertainties in w at the 20% level. 
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3. Cold and hot dark matter 

We measure the density parameter for dark matter to 
be ojd = 0.12 ± 0.01 fairly robustly to theoretical as- 
sumptions, which corresponds to a physical density of 
2.3 x lCT 27 kg/m 3 ±10%. Given the WMAP information, 
SDSS shows that no more than about 12% of this dark 
matter can be due to massive neutrinos, giving a 95% 
upper limit to the sum of the neutrino masses M v < 1.7 
eV. Barring sterile neutrinos, this means that no neutrino 
mass can exceed M v /3 = 0.6 eV. [6] quotes a tighter 
limit by assuming a strong prior on galaxy bias b. We 
show that the recent claim of a neutrino mass detection 
M v > 0.6 eV by Allen et al. hinges crucially on a partic- 
ular low galaxy cluster as measurement and goes away 
completely when expanding the cluster as uncertainty to 
reflect the spread in the literature. 



4- Reionization and astronomy parameters 

We confirm the WMAP team [6] measurement of early 
reionization, t = 0.12+Q;gg. This hinges crucially on 
the WMAP polarization data; using only the unpolar- 
ized WMAP power spectrum, our analysis prefers t = 
and gives an upper limit r < 0.23 (95%). 

Assuming the vanilla model, our Hubble parameter 
measurement h ~ 0.70+q;q 3 agrees well with the HST key 
project measurement h = 0.72±0.07 [49]. It is marginally 
lower than the WMAP team value h w 0.73 ± 0.03 be- 
cause the SDSS power spectrum has a slightly bluer slope 
than that of the 2dFGRS, favoring slightly higher re- 
values (we obtain Q m = 0.30 ± 0.04; the WMAP team 
quote n m = 0.26 ±0.05 [6]). 



E. What have we not learned? 

The cosmology community has now established the ex- 
istence of dark matter, dark energy and near-scale invari- 
ant seed fluctuations. Yet we do not know why they exist 
or the physics responsible for generating them. Indeed, 
it is striking that standard model physics fails to explain 
any of the four ingredients of the cosmic matter budget: 
it gives too small CP-violation to explain baryogenesis, 
does not produce dark matter particles, does not produce 
dark energy at the observed level and fails to explain the 
small yet non-zero neutrino masses. 

Fortunately, upcoming measurements will provide 
much needed guidance for tackling these issues: con- 
straining dark matter properties (temperature, viscosity, 
interactions, etc.), dark energy properties (density evo- 
lution, clustering), neutrino properties (with galaxy and 
cmb lensing potentially sensitivity down to the experi- 
mental mass limits ~ 0.05 eV [137-139]) and seed fluc- 
tuation properties (model-independent measurements of 
their power spectrum [127]). 



The Sloan Digital Sky Survey should be able to make 
important contributions to many of these questions. Red- 
shifts have now been measured for about 350,000 main- 
sample galaxies and 35,000 luminous red galaxies, which 
will allow substantially tighter constraints on even larger 
scales where nonlinearities are less important, as will 
analysis of three-dimensional clustering using photomet- 
ric redshifts [140] with orders of magnitude more galaxies. 
There is also a wealth of cosmological information to be 
extracted from analysis of higher moments of galaxy clus- 
tering, cluster abundance[141], quasar clustering, small- 
scale galaxy clustering[142], Lya forest clustering, dark 
matter halo properties [143], etc., and using information 
this to bolster our understanding the gastrophysics of 
biasing and nonlinear redshift distortions will greatly re- 
duce systematic uncertainties associated with galaxy sur- 
veys. In other words, this paper should be viewed not as 
the final word on SDSS precision cosmology, merely as a 
promising beginning. 
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APPENDIX A: COMPUTATIONAL ISSUES 

In this Appendix, we briefly summarize the technical 
details of how our analysis was carried out. 



1. Monte Carlo Markov Chain summary 

The Monte Carlo Markov chain (MCMC) method is 
a well-established technique [27-29] for constraining pa- 
rameters from observed data, especially suited for the 
case when the parameter space has a high dimensional- 
ity. It was recently introduced to the cosmology commu- 
nity by [31] and detailed discussions of its cosmological 
applications can be found in [7, 32, 33]. 

The basic problem is that we have a vector of cosmo- 
logical data d from which we wish to measure a vector 
of cosmological parameters p. For instance, d might be 
the 1367-dimensional vector consisting of the 899 WMAP 
measurements of the temperature power spectrum Ct for 
I = 2,. ..,900, the 449 WMAP cross-polarization mea- 
surements and the 19 SDSS P(/c)-measurements we use. 
The cosmological parameter vector p might contain the 
parameters of equation (1) or some subset thereof. The- 
ory p is connected to data d by the so-called likelihood 
function £(p,d), which gives the probability distribu- 
tion for observing different d given a theoretical model 
p. In Bayesian analysis, one inserts the actual observed 
data and reinterprets £(p, d) as an unnormalized prob- 
ability distribution over the cosmological parameters p, 
optionally after multiplication by a probability distribu- 
tion reflecting prior information. To place constraints 
on an single parameter, say p-j, one needs to marginalize 
(integrate) over all the others. 

Two different solutions have been successfully applied 
to this problem. One is the grid approach (e.g., [144 
146]), evaluating C(p, d) on a grid in the multidimen- 
sional parameter space and then marginalizing. The 
drawback of this approach is that the number of grid 
points grows exponentially with the number of parame- 
ters, which has in practice limited this method to about 
10 parameters [22]. The other is the MCMC approach, 
where a large set of points Pi, i = l,...n, a chain, is 
generated by a stochastic procedure such that the points 
have the probability distribution C(p, d). Marginaliza- 
tion now becomes trivial: to read off the constraints on 
say the seventh parameter, one simply plots a histogram 
of p 7 . 

The basic MCMC algorithm is extremely simple, re- 
quiring only about ten lines of computer code. 

1. Given p i; generate a new trial point p* = Pi + Ap 
where the jump Ap is drawn from a jump proba- 
bility distribution /(Ap). 



2. Accept the jump (set Pi+i = p*) or reject 
the jump (set Pi+i = Pi) according to the 
Metropolis-Hastings rule [27, 28]: always accept 
jumps to higher likelihoods, i.e., if £(p*,d) > 
£(p i; d), otherwise accept only with probability 
£(p*,d)/£(pi,d). 

The algorithm is therefore completely specified by two 
entities: the jump function /(Ap) and the likelihood 
function £(p*,d). We describe how we compute / and 
C below in sections A 2 and A 3, respectively. 

Table 8 lists the chains we used and their basic proper- 
ties: dimensionality of the parameter space, parameters 
used, data d used in likelihood function, number of steps 
n (i.e., the length of the chain), the success rate (fraction 
of attempted jumps that were accepted according to the 
above-mentioned Metropolis-Hastings rule), the correla- 
tion length (explained below) and the effective length. 
We typically ran a test chain with about 10000 points to 
optimize our choice of jump function / as described in 
Section A 2, then used this jump function to run about 
40 independent chains with different randomly generated 
starting points pi. In total, this used about 30 CPU- 
years of Linux workstation time. Each chain has a pe- 
riod of "burn-in" in the beginning, before it converges to 
the allowed region of parameter space: we computed the 
median likelihood of all chains combined, then defined 
the end of the burn-in for a given chain as the first step 
where its likelihood exceeded this value. Most chains 
burned in within 100 steps, but a small fraction of them 
failed to burn in at all and were discarded, having started 
in a remote and unphysical part of parameter space and 
become stuck in a local likelihood maximum. After dis- 
carding the burn-in, we merged these independent chains 
to produce those listed in Table 8. This standard proce- 
dure of concatenating independent chains preserves their 
Markov character, since they are completely uncorrelated 
with one another. 



2. The jump function / 

As illustrated in Figure 15, consecutive points Pi, i = 
1, ... of a MCMC are correlated. We quantify this by the 
dimensionless autocorrelation function c, shown for the 
reionization parameter r in Figure 16 and defined by 

_ (TjTj+j) - (n) 2 

where averages are over the whole chain. The correlation 
is by definition unity at zero lag, and we define the corre- 
lation length as the number of steps over which the corre- 
lation drops to 0.5. The figure of merit for a chain is its 
effective length N , defined as the number of steps divided 
by the correlation length. Since N is roughly speaking 
the number of independent points, the MCMC technique 
measures statistical quantities such as the standard de- 
viation a p and the mean (p) for cosmological parameters 
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Table 8: Monte Carlo Markov chains used in the chain. The figure of merit for a chain is the effective length (the actual length divided 
by the correlation length). Here we have chosen to tabulate correlation lengths for the T-parameter, since it is typically the largest 
(together with that for n a and u) b , because of the banana degeneracy of Section III A) . The success rate is the percentage of steps 
accepted. "Vanilla" denotes the six parameters (r, uib, u>d, ^A, A s , n s ). In the data column, T denotes the unpolarized power spectrum, X 
denotes the temperaturc/E-polarization cross power spectrum, and r denotes the prior r < 0.3. 



Chain Dim. Parameters 




Data 


Length 


Success Corr. 


length Eff. 


length 


1 


9 Vanilla+f2fc 


+ r + a 


WMAP T+X 


189202 


22% 


218 


868 


2 


7 Vanilla+/^ 




WMAP T+X 


133361 


8% 


78 


1710 


3 


7 Vanilla+ui 




WMAP T+X 


352139 


3% 


135 


2608 


4 


7 Vanilla+n^ 




WMAP T+X 


101922 


7% 


213 


479 


5 


7 Vanilla+r 




WMAP T+X 


178670 


13% 


29 


6161 


6 


6 Vanilla 




WMAP T+X 


311391 


16% 


45 


6920 


7 


6 Vanilla 




WMAP T 


298001 


15% 


25 


11920 


8 


5 Vanilla— n s 




WMAP T+X 


298001 


29% 


7 


42572 


9 


10 Vanilla+H/j 


+ r + a + 


b WMAP T+X + SDSS 


298001 


4% 


69 


4319 


10 


8 Vanilla+/„ 


+ 6 


WMAP T+X + SDSS 


46808 


18% 


24 


1950 


11 


8 Vanilla+iu - 


-b 


WMAP T+X + SDSS 


298002 


4% 


98 


3041 


12 


8 Vanilla+Qjt 


+ 6 


WMAP T+X + SDSS 


298001 


6% 


83 


3590 


13 


8 Vanilla+r + 


6 


WMAP T+X + SDSS 


298001 


12% 


31 


9613 


14 


7 Vanilla+b 




WMAP T+X + SDSS 


298001 


16% 


18 


16556 


15 


7 Vanilla+b 




SDSS+WMAP T 


298001 


16% 


17 


17529 


16 


6 Vanilla — n s 


+ b 


WMAP T+X + SDSS 


298001 


25% 


8 


37250 


17 


6 Vanilla— n s 


+ b 


WMAP T+X + SDSS +0 2 


298001 


25% 


8 


37250 


18 


8 Vanilla+ui - 


-b 


WMAP T+X + SDSS + SN la 


298001 


12% 


25 


11920 


19 


8 Vanilla+r 4 


b 


WMAP T+X + SDSS + SN la 


298001 


5% 


89 


3348 


20 


8 Vanilla+r + 


b 


WMAP T+X + SDSS + r 


151045 


6% 


26 


5809 


21 


8 Vanilla+r + 


b 


WMAP T+X + SDSS + SN la + r 


68590 


6% 


30 


2286 


22 


7 Vanilla+ft 




Other CMB + SDSS 


315875 


30% 


24 


13161 


23 


7 Vanilla+6 




WMAP + other CMB + SDSS 


559330 


20% 


10 


55933 


24 


2 + As 




SDSS 


48001 


41% 


6 


8000 


25 


2 + A s 




SDSS fc max = 0.15 


48001 


36% 


6 


8000 


26 


2 f^m + A s 




SDSS fc max = 0.10 


48001 


31% 


9 


5333 


27 


2 + As 




SDSS no bias corr. 


48001 


38% 


7 


6857 


28 


2 + As 




SDSS linear P(k). 


48001 


50% 


5 


19600 


29 


2 £lm + As 




2dFGRS 


48001 


33% 


9 


5333 




0.2 - 

^ 

1000 2000 3000 4000 

Step 

FIG. 15: The reionization parameter Z as a function of the MCMC 
step. This example is for chain 6 from Tabic 8. 



to an accuracy of order Up/N 1 / 2 . Unless N 1, the 
results are useless and misleading, a problem referred to 
as insufficient mixing in the MCMC literature [29] . 

We attempt to minimize the correlation length by tai- 
loring the jump function to the structure of the likeli- 
hood function. Consider first a toy model with a one- 
dimensional parameter space and a Gaussian likelihood 
C(p) oc e~ p I 2 and a Gaussian jump function f(Ap) oc 
g-Ap /2<r _ what is the best choice of the characteris- 
tic jump size cr? In the limit a — ► oo, all jumps will 
fail; pi = P2 = Cj = 1 for all j and the correlation 
length becomes infinite. In the opposite limit a —> 0, 
almost all steps succeed and we obtain Brownian mo- 
tion with the rms value \pi\ ~ eri 1 / 2 , so it takes of order 
a 2 — > oo steps to wander from one half of the distribu- 
tion to the other, again giving infinite correlation length. 
This implies that there must be an optimal jump size 
between these extremes, and numerical experimentation 
shows that a ~ I minimizes the correlation length. 

In the multiparameter case, strong degeneracies can 
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Steps 



FIG. 16: The autocorrelation function (solid curve) for the exam- 
ple in Figure 15 is seen to be approximately fit by an exponential 
(dashed curve), dropping to 50% at a correlation length of 45 steps 
as indicated by the dotted lines. 



cause a huge correlation length if the jump size is chosen 
independently for each parameter, with the chain tak- 
ing a very long time to wander from one end of the ba- 
nana to the other. A clever choice of parameters that 
reduces degeneracies therefore reduces the correlation 
length. For this reason, the WMAP team used the pa- 
rameters suggested by [26] , and we do the same with the 
minor improvement of replacing A s by A p as the scalar 
normalization parameter as in [24] 5 . To minimize the 
remaining degeneracies, we compute the parameter co- 
variance matrix C = (pp*) — (p)(p)* from the chain it- 
self, diagonalize it as C = RAR', RR* = R*R = I, 
Ay = SijXf, and work with the transformed parame- 
ter vector p' = A _1 / 2 R*[p — (p)] which has the be- 
nign properties (p') = 0, (p'p'*) = I. Inspired by the 
above-mentioned one-dimensional example, we then use 
the simple jump function /(Ap') oc e~' Ap ' / 2cr . We use 
(7 = 1 for all chains except number 1 in Table 8, where we 
obtain a shorter correlation length using a — 0.7. When 
running our test chains to optimize /, we start by guess- 
ing a diagonal C and after the burn-in, we update our 
estimates of both C and the eigenbasis every 100 steps. 
A very similar approach is used in other recent MCMC 
codes, e.g., [32, 147]. 

The WMAP team perform extensive testing to con- 
firm that their chains are properly mixed [7], and we 
have followed the WMAP team in using the Gelman and 



5 When imposing a flatness prior = 0, we retained © s as a free 
parameter and dropped Q^. When additionally imposing a prior 
on h (for the last 6 chains in Table 8), we dropped both © s and 
as f ree parameters, setting = 1 — (w^ + u>i,)/h 2 . 



Rubin i?-statistic [30] to verify that our chains are suf- 
ficiently converged and mixed. Indeed, we find that the 
above-mentioned eigenbasis technique helps further im- 
prove the mixing by cutting our correlation length by 
about an order of magnitude relative to that obtained 
with the WMAP jump function, hence greatly increasing 
the effective length of our chains. 



3. The likelihood function C 

For a detailed discussion of how to compute cosmolog- 
ical likelihood functions, see [7, 148]. Our calculation of 
£(p, d) does little more than combine public software de- 
scribed in other papers, so the details in this brief section 
are merely of interest for the reader interested in exactly 
reproducing our results. 

The total likelihood C is simply a product of likeli- 
hoods corresponding to the data sets used, e.g., WMAP, 
SDSS and SN la. For the CMB, we compute theoret- 
ical power spectra using version 4.3 of CMBfast [23], 
with both the "RECfast" and "PRECISION" options 
turned on. We compute the WMAP likelihood corre- 
sponding to these spectra using the public software pro- 
vided by the WMAP team [7]. Since this software is 
designed for physically reasonable models, not for crazy 
models that may occur during our burn-in, we augment 
it to produce large negative likelihoods for unphysical 
models where it would otherwise give negative x 2_v alues. 
For some of the WMAP+SDSS chains, we evaluate the 
WMAP likelihood Avmap by fitting a quartic polyno- 
mial to ln£ WM Ap from the corresponding WMAP-only 
chains. For this fit, we replace Ud by uj m , u>b by H2, n s 
by H3 and A p by A* inspired by the normal parame- 
ter method of [24]. This approach, described in detail 
in [149], is merely a numerical tool for accelerating the 
computations, and we verify that it has negligible impact 
on our results. 

When combining non-WMAP CMB data with WMAP, 
we include the latest band-power detections from 
Boomerang [150] (madcap), DASI [151], MAXIMA [152], 
VSA [153], CBI (mosaic, even binning) [154] and ACBAR 
[155] with probing effective multipoles I > 600 (where 
they are collectively more sensitive than WMAP) and £ < 
2000 (to avoid complications related to reported small- 
scale excess, which may be due to secondary anisotropics 
or non-CMB effects), which corresponds to the 9+3+3+ 
4 + 6 + 9 = 34 data points plotted in Figure 17. The pre- 
WMAP data has been shown to be consistent both in- 
ternally and with WMAP [156]. We marginalize over the 
quoted calibration uncertainties of 10% for Boomerang, 
4% for MAXIMA and DASI, 5% for CBI, 3.5% for VSA 
and 10% for ACBAR as well as over quoted beam uncer- 
tainties of 15% for Boomerang, 5% for ACBAR and 14% 
for MAXIMA (this last number provides a good fit to the 
combined beam and pointing uncertainties for the three 
measurements used from Table 1 of [152]). We make the 
approximation that all experiments are uncorrelated with 
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FIG. 17: CMB data used. Error bars show do not include the 
calibration and beam uncertainties that we include as described in 
the text. Solid curve corresponds to the "vanilla lite" model of 
equation (2). 



each other and with WMAP, which should be quite accu- 
rate both since sample variance correlations are negligible 
(given their small sky coverage relative to WMAP) and 
since the wmap errors are dominated by detector noise for 
i > 600. When using non-CMB data without WMAP, we 
use all 151 pre- WMAP band power measurements com- 
piled in [157]. For the non-WMAP data, we computed 
the CMB power spectra with the DASh package [158]. 

For SDSS, we compute the likelihood by fitting b 2 times 
the nonlinear power spectrum P(k) to the first 19 band 
power measurements (for k < /c max ~ 0.2/i/Mpc) using 
the window functions and likelihood software provided by 
that paper [20]. For more details about the SDSS data, 
see [159-166]. We compute the nonlinear P(k) using the 
method and software provided by [34]. This software 
takes the linear power spectrum P{k) as input, and we 
compute it using the fitting software provided by [167] for 
the transfer function, the approximation of [39] for the 
linear growth factor and the approximation that P(k) is 
as a product of these two quantities as per equation (C3) 
from [47]. For k > 0.2, this typically agrees with CMB- 
fast 4.3 to better than a few percent. In the absence of 
massive neutrinos (/„ = 0), the separability approxima- 



tion becomes exact and the transfer function fits [167] 
become identical to those of Eisenstein & Hu [168]. 

For SN la, we use the 172 redshifts and corrected peak 
magnitudes compiled and uniformly analyzed by Tonry et 
al. [35] and compute the likelihood with software kindly 
provided by John Tonry. This likelihood depends only on 
(Qa, £^fc, w), and is marginalized over the corrected SN la 
"standard candle" absolute magnitude. 

4. Confidence limits and likelihood plots 

All confidence limits quoted in the tables and text of 
this paper are quantiles, as illustrated in Figure 18. For 
instance, our statement in Table 3 that M v < 1.74eV 
at 95% confidence simply means that 95% of the Mo- 
vables in that chain are smaller than 1.74eV. Simi- 
larly, the entry n s — 0.972^qq 27 in the same column 
means that the distribution of n s -values has median 
0.972, that erfc(2" 1 / 2 )/2 « 15.87% of the values lie 
below 0.972 - 0.027 and that 15.87% of the values lie 
above 0.972 + 0.041, so that 68.27% lie in the range 
n s = 0.9721^027. There is thus no assumption about 
the distributions being Gaussian. In a handful of cases 
involving r, f v and M v , the distribution (see Figure 2 for 
examples) peaks at zero rather than near the median; in 
such cases, we simply quote an upper limit. 

When plotting 1-dimensional distributions f(p) in Fig- 
ure 2, we fit each histogram to a smooth function of the 
form /(p) = e p ^ p ' where P(p) is the 6th order polynomial 
that maximizes the likelihood ]X=i /(P») that the points 
pi,...,p n in the chain are drawn from the distribution 
f(p), subject to the constraint that J f(p)dp = 1. Wc 
found that these smooth curves visually match the raw 
histograms extremely well (see Figure 18 for a rather non- 
Gaussian example) and have the advantage of avoiding 
both the Poisson jaggedness and the excessive smoothing 
inherent in a histogram. 

Our 2-dimcnsional contours are plotted where the 
point density has dropped by e~ Ax I 2 from its maximum, 
where A\ 2 — 6.18 as recommended in §15.6 of [169]. 
These contours would enclose 95% of the points if the 
distribution were Gaussian. When computing the point 
density, there is tradeoff between insufficient smoothing 
(giving contours dominated by Poisson noise) and exces- 
sive smoothing (which artificially broadens the contours, 
particularly in the narrow direction of a degeneracy ba- 
nana). We found that this was alleviated by computing 
the contours in the linearly transformed 2-dimensional 
space defined in Section A 2 where the covariance matrix 
is the identity matrix. 
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